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ABSTRACT 

We report on new VUT optical spectroscopic and multi- wavelength archival observations of SN 1996cr, a 
previously identified ultraluminous X-ray source known as Circinus Galaxy X-2. The spectrum of the optical 
counterpart confirms SN 1996cr as a bona fide type Iln SN, in accordance with its tentative SN classification at 
X-ray wavelengths. SN 1996cr is one of the closest SNe (w3.8 Mpc) in the last several decades and in terms of 
flux ranks among the brightest radio and X-ray SNe ever detected. Optical imaging from the Anglo- Australian 
Telescope archive allows us to isolate the explosion date to between 1995-02-28 and 1996-03-16, while the 
wealth of optical. X-ray, and radio observations that exist for this source provide relatively detailed constraints 
on its post-explosion expansion and progenitor history, including an preliminary angular size constraint from 
VLBI. Archival X-ray and radio data imply that the progenitor of SN 1996cr evacuated a large cavity just prior 
to exploding via either a sped-up wind or a pre-SN explosion. The blast wave likely spent ^1-2 yrs in relatively 
uninhibited expansion before eventually striking the dense circumstellar material which surrounds SN 1996cr 
to become a prodigious X-ray and radio emitter The X-ray and radio emission, which trace the progenitor 
mass-loss rate, have respectively risen by a factor of >2 and remained roughly constant over the past « 7 yr. 
This behavior is reminiscent of the late rise of SN 1987 A, but three orders of magnitude more luminous and 
much more rapid to onset. SN 1996cr may likewise provide us with a younger example of SN 1978K and 
SN 1979C, both of which exhibit flat X-ray evolution at late times. The optical spectrum suggests that the 
progenitor was a massive star that shed some of its outer envelope (many Mq) prior to explosion, while the 
complex Oxygen line emission hints at a possible concentric shell or ring-like structure. Taken together, this 
implies that a substantial fraction of the closest SNe observed in the last several decades have occurred in 
wind-blown bubbles, and argues for the phenomena being widespread. 

Subject headings: stars: supernovae: general — stars: circumstellar matter — X-rays: supernovae — radio: 
supernovae — 



1. INTRODUCTION 

Although thousands of supernovae (SNe) have been discov- 
ered to date, only several dozen have been detected at either 
X-ray^ or radio"^ wavelengths (designated XSNe or RSNe, re- 
spectively). Many of these XSNe and RSNe were originally 
detected simply as variable X-ray or radio sources, and only 
verified as true SNe after careful examination of archival opti- 
cal data or through optical follow-up (e.g., Rupen et al. 1987; 
Ryder etal. 1993). This particular path to discovery has meant 
that only a handful of these sources have been well-studied 
during the first several hundred days, a period which is critical 
for identifying and characterizing the true nature of the SNe. 
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By contrast, XSNe and RSNe do comprise a substantial fraction of all 
SNe well-studied beyond a few hundred days, as such SNe end up being 
nearby and remain relatively bright for years to decades. 



Such SNe typically turn out to be core-collapse SNe of ei- 
ther type Ibc, which have been associated with long-duration 
gamma-ray bursts (GRBs), or type II, whereby the intense 
X-ray and/or radio emission is thought to stem from the inter- 
action between the expanding shock and a dense progenitor 
wind. When well-sampled, the additional multi-wavelength 
constraints from XSNe and RSNe can provide physical in- 
sights into the late evolutionary stages of massive stars that 
are otherwise impossible to obtain. For instance, while ro- 
bust constraints on the overall time-averaged mass-loss in var- 
ious phases of stellar evolution exist, there has been a long- 
standing debate over the number and sequence of various evo- 
lutionary stages for massive stars, and relatively few observa- 
tional constraints on the actual evolution of mass-loss within 
many particular stages (e.g., Lamers et al. 1991; Langer et al. 
1994; Stothers & Chin 1996; Maeder & Meynet 2000). Be- 
cause the SN blast wave travels orders of magnitude faster 
than the stellar wind, studying the interaction between the 
blast wave and the progenitor stellar wind allows us to probe 
tens of thousands of years of evolution in a matter of decades. 
SNe progenitors are considered to provide the bulk of all pro- 
cessed stellar material additionally, and thus the characteriza- 
tion of these sources can likewise help further our understand- 
ing of overall galactic chemical enrichment. 

Here we report on the spectroscopic confirmation of one 
such serendipitous source, SN 1996cr (Bauer 2007), and 
present a multi-wavelength follow-up study using archival 
data to determine its explosion date and temporal properties. 
SN 1996cr was originally detected as Circinus Galaxy (CG) 
X-2, an ultraluminous X-ray source in the nearby Circinus 
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Galaxy (Freeman et al. 1977), which Bauer et al. (2001) found 
to exhibit many characteristics of a young, type II SN en- 
shrouded in a dense circumstellar environment. In particu- 
lar, it demonstrated a factor of >30 increase in X-ray flux 
between 1997-2000, a fcT ~10 keV thermal spectrum with a 
strong, blended Fe emission-Une component at 6.85 keV, and 
spatial coincidence with a strong radio and Ha-emitting point 
source. The proximity of this SN affords us a rare opportunity 
to study in detail a type II SN which is strongly interacting 
with its circumsteUar medium. 

This paper is organized as follows: data and reduction 
methods are detailed in §2; confirmation of SN 1996cr as a SN 
and isolation of its explosion date are provided in §3; overall 
temporal and spectroscopic constraints for SN 1996cr are in- 
vestigated in §4; and finally conclusions and future prospects 
are explored in §5. Throughout this paper, we adopt a dis- 
tance of 3 .8 ± 0.8 Mpc to the Circinus Galaxy (converted from 
Freeman et al. 1977). While the Circinus Galaxy lies close to 
the Galactic Plane {\b\ =3.°8), it is located within a Galactic 
"window" with a visual extinction of Ay = 1 .5 ± 0.2 and a 
neutral hydrogen column density of A^h = (3.0 ±0.3) x 10^' 
cm~^ (whereas neighboring regions typically have Ay = 3.0 
and A^H = (5-10) x lO^' cm'^; Schlegel et al. 1998; Dickey & 
Lockman 1990). Due to its ~ 65° inclination, however, there 
is significant internal obscuration as well [A^h ~ (5-8) x 10^' 
cm~^ typically; Bauer et al. 2001]. 

2. OBSERVATIONAL DATA AND REDUCTION METHODS 

We describe below the extensive observational data used 
to constrain the properties of SN 1996cr. Fortunately, plen- 
tiful archival data exist for the Circinus Galaxy due to the 
fact that it hosts the second closest Compton-thick active 
galactic nuclei (AGN) to our own Galaxy and exhibits signs 
of vigorous star formation (e.g., Matt et al. 1999). The 
resulting multiwavelength dataset provides good constraints 
on the explosion date and yields useful long-term X-ray 
and radio light curves. SN 1996cr lies 25"to the south of 
the Circinus Galaxy nucleus. We adopt the position of 
a = 14M3'"10?01 S = -65°20'44:'4 (J2000), determined from 
the radio observations. 

2.1. VLT Spectroscopy 

SN 1996cr was visited several times in service mode with 
the VLT FORS I spectrograph in two separate campaigns. 
The first program was initiated over the course of 2005-03- 
06 through 2005-03-17 but never completed, resulting in a 
low signal-to-noise spectrum. As such, we do not provide any 
further details. The second program was executed over three 
nights (2006-01-26, 2006-02-02, and 2006-03-10), yielding 
a high-quality spectrum confirming the ambiguous features 
seen in the spectrum from our first program. Our analyses 
focus only on the spectrum from the first night, as the obser- 
vations on the subsequent two nights suffer from significantly 
worse seeing («l."5-2"0) such that the co-addition of these 
frames failed to improve the signal-to-noise of the final spec- 
trum. Our observations on the first night consisted of four 
1200 s exposures taken with the 300V grism and the TEK 
(24/x pixel) CCD using a slitwidth of 2f'0. The grism and 
CCD combination provided a dispersion of 2.66 A pixel"' and 
a total useful wavelength coverage of 3600-7990A. These 
images were taken in excellent seeing conditions, in which 
stellar sources had a full-width half-maximum (FWHM) mea- 
sured at 0."6-0"7 on the FORS I chip over the course of the 
night. Our target, SN 1996cr, had an w 20% larger FWHM 



image size than the point sources in the acquisition image, 
suggesting that it may be spatially extended or possibly con- 
taminated by a coincident H ii region. We investigated this 
further by measuring the FWHM orthogonal to the disper- 
sion axis for the strong emission fines [O III]A5007 and Ha, 
which had sizes of 0."74 and 0."85, respectively. The in- 
herent FWHM slowly increased with decreasing wavelength, 
such that the FWHMs quoted above were 5±3% and 35±4% 
larger than their stellar equivalents, respectively. Thus the 
[O III]A5007 FWHM is formally unresolved, while the Ha 
line is extended. The strength, extent, and visible asymmetry 
in the FWHM profile of the extended Ha are fully consistent 
with the underlying H II region (see §3 for details), although 
we caution that marginally-resolved light echoes could po- 
tentially contaminate and broaden the PSF (e.g., Sugerman 
2003). The spectral resolution (11 A, or 500km s"') was es- 
timated from the width of the arclines taken with a 1 ."0 slit, 
which matched closely the on-sky image size of SN 1996cr. 
After bias subtraction and flat-fielding using standard tech- 
niques, the co-added spectra were extracted, wavelength cal- 
ibrated, and flux calibrated (using the standard LTT4816). 
Fig. 1 shows the resulting optical spectrum from 2006-01-26. 
We conservatively estimate the flux errors on the spectrum at 
20%. 

2.2. Archival Optical Data 

We searched through the European Southern Observatory 
(ESO), Anglo Australian Telescope (AAT), United Kingdom 
Schmidt Telescope (UKST), and Hubble Space Telescope 
(HST) archives. We only list here data that ultimately were 
used to constrain the explosion date of SN1996cr. Rele- 
vant data are listed in Tables 1-4. Aperture photometry was 
performed throughout, with PSF-fitting employed as a cross- 
check; there is only one instance where the PSF-fitting mag- 
nitude differed from the aperture-measured value at > la, 
which we explicitly document below. 

2.2.1. ESO 

The Circinus Galaxy was observed extensively with SUSI 
at the NTT on 1993 April 9 using several narrow-band fil- 
ters (#369:[0 III|, #430:5100A continuum, #629:HaH-[N II], 
#700:[S II], #443:7000A continuum, #415:[Fe XI]A7892)" 
as weU as with die IRAC2 camera at flie ESO/MPI 2.2m tele- 
scope on 1994 June 25 using broad-band J and H filters. 
Reduced images were kindly provided by E. Oliva and A. 
Marconi (private communication, 2007), for which cutouts of 
SN 1996cr are shown in Fig. 2. Details of the observations 
and reduction procedures are given in Marconi et al. (1994). 
The [O III], Ha-H[N II], and [S II] images were the only im- 
ages provided to us in a flux-calibrated state, and therefore are 
the only ones for which we measure photometry. Magnitudes 
for the region in the vicinity of SN 1996cr were measured us- 
ing a 1 "0 radius circular aperture and previously established 
zero points. An aperture correction of 0.2 mags was estimated 
empiricaUy using several bright, isolated point sources in the 
images. The resulting aperture-corrected magnitudes are pre- 
sented in Table 1. 

2.2.2. HST 

The Circinus Galaxy was observed with both the WFPC2 
and NICMOS instruments aboard HST on four separate oc- 
casions, as outlined in Table 2. Details of the photometric 

' ' http://filters.ls.eso.org/efs/efs_fi.htm 
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Fig. 1. — Flux-calibrated VLT FORS spectrum of SN 1996cr from 2006 January 26. The spectrum suffers from an extinction of Ay ~ 4-6, which we 
have not coirected. Obvious emission lines have been identified. There is strong, naiTow Ha emission (FWHMfa 714 km s"') as well as several broad O 
complexes showing a distinct underlying velocity structure. An empirical model has been fit to the data as described in § 4.2; the green line shows the overall 
fit to the spectrum while the red (blue) lines show the individual nan'ow (broad) components. Strong, broad unidentified residuals remain around the He I lines 
(AA5876,7065); this emission is likely associated with the structure of a He shell or the progenitor wind, but is not easily disentangled. 



and astrometric reduction of the WFPC2 data are given in 
Bauer et al. (2001). Briefly, after standard calibration of the 
HST images, we used the IRAF package DAOPHOT to mea- 
sure aperture-corrected magnitudes using a 0."2 radius aper- 
ture for all sources down to the la detection limit. The HST 
images were then aligned to the Hipparcos/Tycho astromet- 
ric reference frame to rs0"4. Due to the negligible overlap 
of the F606W observation with the other filters and complete 
lack of coverage of SN 1996cr itself, we do not discuss it fur- 
ther For the NICMOS data, of which only the N1C3 images 
provided useful imaging, we used the standard pipeline data 
products. We performed photometry using a 0"5 radius aper- 
ture and applied band-dependent aperture corrections deter- 
mined from a set of isolated, unsaturated point sources in the 
vicinity of SN 1996cr. Cutout images of SN 1996cr are pre- 
sented in Fig. 3. Notably, SN 1996cr appears to lie at the cen- 
ter of a diffuse patch of Ha emission in the F656N image 
and our aperture magnitude here differs from a PSF-fitting 
one at > 3cr (0.11 magnitudes). We thus adopt the PSF-fitting 
magnitude, as well as an additional systematic error of 0.1 
magnitudes to reflect the larger uncertainty associated with 
the deblending. 

2.2.3. AAT 

The Circinus Galaxy was observed numerous times with the 
TAURUS Fabry-Perot instrument on the AAT during 1995- 
02-21 to 1995-02-28 and 1996-03-15 to 1996-03-20. The 



data were retrieved through the AAT archive.'^ The ob- 
servations were comprised primarily of narrow-band imag- 
ing and TAURUS spectral imaging cubes centered on the 
He IIA4686, [O III]AA5007,4959, Ha+[N II]AA6583,6548, 
and [S II]AA6731,6716 emission lines. The original investi- 
gators discuss their observational design in Veilleux & Bland- 
Hawthorn (1997) and Elmouttie et al. (1998), while more 
complete data reduction procedures are outlined in Gordon 
et al. (2000). To summarize, the TAURUS instrument was 
used in the angstrom imaging mode, wherein narrow-band 
filters were used at different tilt angles to isolate the lines of 
He II, [O III], Ha+[N II], and [S II] (Bland-Hawthorn & Jones 
1998). Each square pixel subtended 0."315 on the sky and the 
atmospheric seeing at FWHM averaged «l"2-2."0. Unfor- 
tunately, a significant portion of the archived TAURUS im- 
ages for the Circinus Galaxy have only limited accompanying 
CCD calibration data; there were often no obvious superbias 
frames, dark frames, dome-flats, or flux standards, so the pho- 
tometric quality of the reduced data is limited. 

Our best-effort reduction proceeded as follows. The over- 
scan regions were used to subtract the bias from each frame 
and each image was flat-fielded using a sky flat. When mul- 
tiple frames in a given filter were available, the images were 
combined to reject cosmic rays using the IRAF task CRREJ. 
Source detection and photometry were performed with SEx- 
tractor (Bertin & Arnouts 1996). A set of 40 relatively bright, 

http://site.aao.gov.au/arc-bin/wdb/aat_database/observation_log/make 
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Fig. 2. — Postage-stamp images centered on SN 1996cr from ESO imaging observations in 1993-1994 courtesy of E. Oliva and A. Marconi (private commu- 
nication, 2007). Images are 24" on a side. The position of SN 1996cr is denoted by a 2" radius circle. Each image is labeled, denoting both the instrument and 
filter used to acquire it and the observation date. 
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Fig. 3. — Postage-stamp images centered on SN 1996cr from //ST WFPC2 and N1C3 imaging observations spanning 1996 to 1999. SN 1996cr appears as a 
strong, unresolved emission-line source in the F502N and F656N images, but displays only weak continuum as evidenced its unremarkable appearance in all of 
the other //Sr images and the fainter broad-band magnitudes in Table 2. See Fig. 2 caption for details of image properties. 




Fig. 4. — Postage-stamp images centered on SN 1996cr from AAT imaging observations. The top six images from 1995-02-21 to 1995-02-28 represent the 
comparison epoch, while the bottom six images from 1996-03-16 to 1996-03-19 represent the discovery epoch. Despite the poor quality of some images, it is 
clear that SN 1996cr is consistently luminous in all discovery images spanning several days. See Fig. 2 caption for details of image properties. 
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Fig. 5 . — Postage-stamp images centered on SN 1996cr from UKST imag- 
ing observations in 1998. SN 1996cr is still evident as a slight excess in the 
Ha images compared to similar narrow-band images in Figs. 2 and 4. See 
Fig. 2 caption for details of image properties. 

isolated, unsaturated stars, which are detected in all the AAT 
and HST images, were chosen to register the AAT images to 
the HST coordinate frame. The aUgnment of each AAT im- 
age to this reference frame is accurate to ?s0."l-0."2 (la). 
The zeropoints for the narrow-band images were obtained 
by bootstrapping the AAT magnitudes to the well-determined 
HST ones. To this end, HST colors were obtained for the 40 
stars and matched to standard main sequence stellar templates 
from the Bruzual Atlas, modified by a standard Galactic 
dust model with E(B - V) = 1 (the extinction measured to- 
ward the Circinus Galaxy; e.g., Schlegel et al. 1998) using 
the IRAF package SYNPHOT. Among the 40 calibration stars, 
25 had colors consistent with one of the stellar templates. 
These particular stellar templates and the TAURUS-specific 
transmission curves'^ were input into SYNPHOT to convert be- 
tween HST and TAURUS magnitudes. We used the iRAF task 
FlTPARAMS to compare interactively the i/5r-derived magni- 
tudes to the SExtractor ones and derive reasonable zeropoints 
for each image. We typically rejected a few outUers in fit- 
PARAMS to obtain an adequate fit (0.3-0.5 mag Icr error typ- 
ically on the zeropoint). The He II, [O III], and [S II] bands 
were observed multiple times, so we combined zeropoint esti- 
mates to ensure that magnitudes for the detected sources were 
consistent to within «0. 1-0.2 mags across these images. For 
bands which sampled comparable spectral windows, we used 
SYNPHOT to convert from one band to another (for instance, 
from TAURUS 5020 A/30 A to HST ¥502^ or from the VLT 
FORS spectrum to all overlapping bands). The resulting mag- 
nitudes are listed in Table 3, while Fig. 4 shows 12 images 
taken in 1995 February and 1996 March. 

2.2.4. UKST 

The Circinus Galaxy was observed on numerous occasions 
with the UKST as detailed in Table 4. The data were retrieved 
from the SUPERCOSMOS archive, although only a subset 
of the photographic plates have been properly digitized and 
archived in a usable form. Fig. 5 shows four images taken in 
1998 March and April. The Ha images from 1998 (as well 

http://www.stsci.edu/hst/observatory/cdbs/bz77.htmI 
http://www.aao.gov.au/local/www/cgt/ccdimguide/filtercat.html 
" http://surveys.roe.ac.uk/ssa/ 



as similar degraded ones from 1999 and 2001) show a slight 
enhancement over the 1995 Ha reference images from the 
AAT and ESO. For the four images that were taken after the 
discovery window of SN 1996cr and properly digitized, we 
performed photometry using a If'O radius circular aperture 
and previously established zero points. An aperture correction 
of 0.2 mags was estimated empirically using several bright, 
isolated point sources in the images. The resulting aperture- 
corrected magnitudes are presented in Table 4. 

2.2.5. Swift UVOT 

The Circinus Galaxy was observed on 2007-03-23 with the 
Swift Ultraviolet/Optical Telescope (UVOT) for 617-2461 s 
depending on the filter. The images were reduced using stan- 
dard pipeline procedures and photometry was performed on 
the images using a 3."0 radius circular aperture. SN 1996cr 
was not detected in any of the six UVOT filters, with aperture- 
corrected 3cr upper hmits of V > 19.0, B > 20.2, U > 20.0, 
UVWl > 20.4, UVW2 > 20.8, and UVM2 > 21.2. 

2.2.6. SNe Optical Monitoring 

There are three additional useful upper-limit constraints 
based on the SNe monitoring observations of the Circinus 
Galaxy by Rev. Robert Evans on 1995-03-31, 1995-08-13, 
and 1996-03-03 (private communication, 2000). We have 
adopted upper limits of V > 14 based on R. Evans' estimated 
limiting magnitude ofV^i 15. 

2.3. Archival X-ray Data 

Details of the various X-ray observations are given in Ta- 
ble 5. We describe below our reduction methods for each 
dataset. When possible, the X-ray fluxes and absorption- 
corrected luminosities for SN 1996cr were calculated via 
spectral analysis in XSPEC (Arnaud 1996) using the Cash 
statistic (Cash 1979). Unless stated otherwise, errors on spec- 
tral parameters are for 68% confidence, assuming one pa- 
rameter of interest. Following the treatment of SN 1987 A 
(e.g.. Park et al. 2005), we characterized the X-ray spec- 
tra of SN 1996cr with an absorbed variable-abundance non- 
equilibrium ionization (NET) shock model (vpshock using 
NEI version 2.0 in XSPEC) based on ATOMDB (Smith et al. 
2001)."' The best-fitted parameters to the combined XMM- 
Newton and Chandra HETGS 0th order dataset were as fol- 
lows: A^H = 8.4 X 10^' cm-2, kT = 13.4 keV, Zsi = 2.25Z0^si, 
Zs = 3.5Zq,s, -Zca = 8ZQ_ca, and Zpe = 2ZQ,Fe. We note, 
however, that these specific abundances, and even the esti- 
mated temperature to some extent, should be used with cau- 
tion, as there are known to be significant deviations between 
adopted models (e.g., Nymark et al. 2006). The typical col- 
umn density, A^h, derived from spectral fits to the data is 
much larger than the estimated Galactic absorption colunm 
(«3 X 10^' cm"^), implying significant internal absorption ei- 
ther from the disk of the Circinus Galaxy at large or the imme- 
diate vicinity of SN 1996cr. Fluxes were measured by vary- 
ing the normalization of the above model in the 0.5-2 keV or 
2-8 keV bands. Our adopted model fit all of the spectral data 
adequately, although we note that there was typically some 
excess residual emission around prominent emission lines, as 
well as an apparent ~ 2cr soft excess in the 2004 HETGS data. 

While this model is believed to best characterize the X-ray emission 
from young SNe such as SN 1996cr, there are still many serious caveats. See 
http://cxc.harvard.edu/atomdb/issues_caveats.html for details. 
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TABLE 1 
ESO Observations 



Instrument 


Filter 


Date 


Exp. 


Seeing 


Magnitude 


SUSI/NTT 


ESO#369 ([O ni]) 


1993-04-09 


900 


o;'7 


>21.7 


SUSI/NTT 


ESO#430(5108Acont.) 


1993-04-09 


900 


o;'7 




SUSI/NTT 


ESO#629 (Ha+[N II]) 


1993-04-09 


480 


o;'7 


19.1±0.2t 


SUSI/NTT 


ESO#700 ([S II]) 


1993-04-09 


900 


o;'7 


21.2±0.3t 


SUSI/NTT 


ESO#443(7027Acont.) 


1993-04-09 


900 


o;'7 




SUSI/NTT 


ESO#415([FeXI]A7892) 


1993-04-09 


900 


o;'7 




IRAC2/ESO-MPI 2.2m 


J (1.25 /im) 


1993-06-25 


360 


o;'9 


>22.0 


IRAC2/ESO-MPI 2.2m 


H (1.65 /im) 


1993-06-25 


360 


o;'9 


> 19.8 



Note. — Column 1: Instrument and telescope. Column 2: Filter. Column 3: UT date of observation given as 
year-month-day. Column 4: Exposure time in seconds. Column 5: Seeing. Column 6: Aperture-corrected Vega 
magnitude or 3o- upper limit. Details are given in § 2.2.1. The detected magnitudes denoted by fs provide constraints 
on line emission trom the H II region which spatially overlaps SN 1996cr. 



TABLE 2 
HST Observations 



TABLE 3 
AAT TAURUS Observations 



Instrument 


Filter 


Date 




Exp. 


Magnitude 


WFPC2 


F606W 


1996-08- 


■11 


200, 400 




NIC3 


F164N 


1998-03- 


16 


64 


20.75 ±0.28 


NIC3 


F166N 


1998-03- 


16 


64 


> 21.09 


NIC3 


F196N 


1998-03- 


16 


80 


21.04 ±0.29 


NIC3 


F200N 


1998-03- 


16 


80 


20.59 ±0.27 


NIC3 


F212N 


1998-10- 


16 


160 


20.90 ±0.1 8 


NIC3 


F215N 


1998-10- 


■16 


160 


20.60 ±0.1 9 


WFPC2 


F502N 


1999-04- 


■10 


900, 900 


20.67 ± 0.06 


WFPC2 


F547M 


1999-04- 


■10 


60 


22.81 ±0.12 


WFPC2 


F656N 


1999-04- 


■10 


800, 800 


17.98±0.11 


WFPC2 


F814W 


1999-04- 


■10 


40 


21.20 ±0.09 



Note. — Column 1: /f5r instrument and filter Filter denotes the central 
filter wavelength in nanometers and the filter width (N-narrow, M-medium, 
W-wide). Details on individual filters can be found in the WFPC2^' and NIC- 
MOS'^ Instrument Handbooks. Column 2: UT date of observation given 
as year-month-day. Column 3: Exposure time(s) in seconds. Column 4: 
Aperture-corrected Vega magnitude or 3o- upper limit. Magnitudes were mea- 
sured with circular apertures of radius 0"2 and 0"5 for WFPC2 and NIC3, 
respectively, and corrected for missing flux due to the shape of the PSF. Note 
that increasing the aperture radius in the HST F656N band to l"2 yields a 
magnitude of 17.5 (i.e., an increase of ^25% over the point-like magnitude 
from SN 1996cr alone) which we attribute to the fiux of the underlying H 11 
region. We have added a systematic error of 0.1 magnitude to the F656N 
value to reflect the accuracy to which we can deblend the point source emis- 
sion from the underlying H II region. The H II region alone should have a 
magnitude of 18.7 ± 0.2, which is consistent with values measured from our 
early-time comparison images. Likewise, increasing the aperture in the NIC3 
observations any further will begin to include nearby contaminating point 
sources, as shown in Fig. 3. Finally, we have incorporated an additional sys- 
tematic error of 0.1 mag to account for the intrapixel sensitivity variations in 
the NIC3 observations. 
°http://www.stsci.edu/instruments/wfpc2AVfpc2_hand_current/ 
*http://www.stsci.edu/hst/nicmos/documents/handbooks/v4.1/ 

These issues hiave a negligible effect, however, on our flux es- 
timates to within errors and will be addressed in more detail in 
a separate pubhcation on X-ray spectral analysis of SN 1996cr 
(F. Bauer et al., in preparation). 

2.3.1. Chandra 

The Circinus Galaxy was observed on several occa- 
sions with Chandra using the backside-illuminated Advanced 
Imaging CCD Spectrometer (ACIS-S) in the focal plane, 
both with and without the High-Energy Transmission Grat- 
ing Spectrometer (HETGS). Archival ACIS-S and HETGS 
zeroth-order data were retrieved from the Chandra Data 
Archive and processed following standard procedures using 
CIAO (v3.4) software. Additionally we removed the 0"5 
pixel randomization, corrected for charge transfer inefficiency 





Fi: 


Iter 


5020/30 


(Of 


=0) 


6583/45 


(Of 


=9) 


6583/45 


(9f 


=9) 


6618/55 


{Of 





Date 



Exp. .Seeing Magnitude 



R 
B 

4690/26 (Of =5) 
6730/30 (Of =5) 
4690/26 Z-CUBE (fif =5) 
4690/26 (Of =0) 
4690/26 (fif =0) 
5020/30 (9f =0) 
6730/30 (Of =0) 



1995-02-21 
1995-02-22 
1995-02-28 
1995-02-28 
1995-02-28 

1995- 02-28 

1996- 03-16 
1996-03-16 
1996-03-16 
1996-03-17 
1996-03-19 
1996-03-19 
1996-03-19 



300 
300 
267 
1442 
60 
20 
1200 
1200 

1200 
1000 
120 
1000 



1"2 
lf'2 

i;'i 
ir'i 

1^2 

i;'4 
3;'5 
3;'o 
3;'o 
3;'i 
i;'2 
i;'3 
i;'3 



> 20.0 

> 18.3 

> 18.1 

> 18.4 

> 19.6 
>22.1 
19.5 ±0.6 
15.8 ±0.5 

19.8 ±0.5 
19.7 ±0.3 
18.7 ±0.4 
15.4 ±0.3 



Note. — Column 1: Filter, given here as central wavelength and bandwidth in 
Angstroms. 6370/30, 6583/45. 5()2()/3(), 4690/26 are narrow-band filters used with 
the AAT TAURUS Tunable Filter (see AAT Filter Catalog" for details), while B and 
R are standard Johnson filters. Note that for the TAURUS instrument, we also list a 
tilt angle, which tunes the wavelength range such that \^ff » A„„„(l -eJ/28984), 
where Of is tilt angle in degrees. Column 2: UT date of observation given as year- 
month-day. Column 3: Exposure time in seconds. Column 4: Seeing. Column 5: 
Aperture-corrected Vega magnitude or 3(T upper limit. Details are given in § 2.2.3. 

" http ://w w w. aao . gov. au/local/w w w/cgt/ ccdimguide/filtercat.html 

(CTI), performed standard ASCA grade selection, excluded 
bad pixels and colunms, and screened for intervals of exces- 
sively high background (none was found). Analysis was per- 
formed on reprocessed Chandra data, primarily using CIAO, 
but also with FTOOLS (v6.3) and custom software including 
ACIS EXTRACT (v3.128; Broos et al. 2005). Spectra were ex- 
tracted using ACIS EXTRACT with a 95% encircled-energy re- 
gion derived from the Chandra PSF library. As SN 1996cr 
lies at the edge of diffuse emission associated with the AGN, 
a background spectrum was extracted from a local annular re- 
gion after excluding nearby point sources and the strongest 
portions of the circumnuclear halo and ionization cone asso- 
ciated with the AGN (the annulus size was grown until ap- 
proximately 100 counts could be extracted). We generated 
calibration products (including an aperture correction) for the 
HETGS/ACIS-S spectra and fit them separately within xspec 
using our best-fit model to estimate the flux. The majority of 
the observations were performed in configurations which mit- 
igate pile-up; when significant we have noted the estimated 
pile-up fractions for SN 1996cr, and corrected for them using 
the XSPEC pileup model of Davis (2001). 

2.3.2. XMM-Newton 
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TABLE 4 
UKST Observations 



Plate # 


Filter 


Date 


Exp. 


Grade 


Magnitude 






1 QQ 1 07 17 


2400 


AI3* 




OR 1 1 1 


±±±£ii t\j\j^y\j 


1 QQ4-0(S-97 


2480 


BI3* 




OR16238 


IIIaF/OG59() 


1994-08-12 


3300 


AID2 




OR17455 


IIIaF/OG59() 


1997-03-11 


300 


al* 




OR 17484 


IIIaF/OG59() 


1997-03-31 


300 


a* 




HA17930 


4415/HA659 


1998-03-05 


10800 


A2 


17.1 ±0.6 


OR17954 


4415/00590 


1998-03-19 


900 


a 


18.3 ±0.9 


HA 17987 


4415/HA659 


1998-04-05 


10800 


AT2 


17.3 ±0.4 


OR17988 


4415/OG590 


1998-04-05 


900 


a 


18.0 ±0.7 


OR18316 


4415/OG590 


1999-03-11 


900 


a* 




HA18323 


4415/HA659 


1999-03-12 


2160 


bU* 




HA18328 


4415/BiA659 


1999-03-14 


10800 


bT* 




HA19197 


4415mA659 


2001-05-24 


8400 


AU2* 





Note. — Column 1: A one- or two-character code for the filter/emulsion combi- 
nation along with a running number for all UKST plates. Column 2: Emulsion and 
filter. OG590 is red, while HA659 is Ha. Column 3: UT date of observation given 
as year-month-day. Column 4: Exposure time in seconds. Column 5: All plates are 
quahty controlled and assigned a grade. The first letter of the grade (usually A,B7C) 
indicates the overall quality of the plate. The subsequent letters indicate specific de- 
fects as follows: I - Denotes image size larger than 40 microns. T - Denotes detectable 
image elongation on most images. U - Denotes underexposure (relative to exposure 
time). D - Denotes overexposed plate (high central density). Survey plates are nor- 
mally also given a grade number; the lower the number the better the plate quality. 
An 'A' grade plate scores 3 or lower Plates denoted by '*' were only available as 
lower quality "Finder" (GIF) format and thus are not as photometrically reliable. Ad- 
ditionally, the spatial resolution of plate OR16238 was too poor to provide a usefiil 
comparison magnitude due to source blending. Column 7: Aperture-corrected Vega 
magnitude. Details are given in § 2.2.4. 

The Circinus Galaxy was observed once withXMM-Newton 
using the EPIC p-n and MOS1/MOS2 detectors on 2001-08- 
06 for « 104 ks. Archival data were retrieved from the XMM- 
Newton Science Archive and processed following standard 
procedures with SAS (v7.0.0). Additionally, we screened 
for and removed intervals of excessively high background (» 
15% for MOS1/MOS2 and « 40% for p-n), leaving 85.5 ks, 
91.8 ks, and 59.5 ks of useful exposure with the MOSl, 
MOS2, and p-n instruments, respectively. Spectra were ex- 
tracted with SAS using a circular aperture of radius 1 1 ."0 and 
local background was extracted in a manner similar to that 
done for the Chandra observations. We generated calibration 
products for the p-n/MOSl/MOS2 datasets and joint-fit the 
three spectra within XSPEC using our best-fit model to esti- 
mate the flux. 

Additionally, we used the XMM-Newton data to place upper 
limits on the potential flux from SN 1996cr in the ASCA and 
BeppoSAX datasets, both of which lack the spatial resolution 
to resolve the galaxy into its various components (e.g., see 
Fig. 1 of Bauer et al. 2001). To this end, spectra for the en- 
tire Circinus Galaxy were extracted from a 2f5 radius circular 
aperture, excluding emission from SN 1996cr. Unfortunately, 
the ultraluminous X-ray source CG X-1 is known to vary with 
a period of 7.5 hr and comprises a non-negligible fraction of 
the total X-ray emission from Circinus. Thus we extracted 
a pessimistic spectrum additionally excluding CG X-1. This 
provides conservative upper limits for the total emission from 
the Circinus Galaxy during the ASCA and BeppoSAX obser- 
vations, under the assumption that the fluxes of the various 
emission components within the Circinus Galaxy (aside from 
CG X-1 and SN 1996cr) have not changed since 1995. Based 
on the spatially resolved observations that we do have, this 
assumption appears reasonable. More optimistic upper lim- 
its can be estimated by further subtracting the time-averaged 
flux from CG X-1. Given ambiguities in determining how the 



unresolved observations sample the light curve of CGX-1, 
we simply provide here the flux of CG X-1 in its "high" state 
(see Bauer et al. 2001) which is 5.7 x 10"'^ ergs s"' cm~-^ 
and 4.3 x lO'^^ ergs s""' cm'^ in the 0.5-2 keV and 2-8 keV 
bands, respectively, and note that optimistic upper limits for 
the flux of SN 1996cr would likely be a factor of «2-3 lower. 

2.3.3. ROSAT 

The Circinus Galaxy was observed on several occasions 
with ROSAT using the High-Resolution Imager (HRI). Pro- 
cessed data were retrieved from the High Energy Astrophysics 
Science Archive Research Center (HEASARC), and analy- 
sis was performed using ftools and custom software. We 
adopted the best-fit model from the combined observations 
(see above) to estimate the ROSAT HRI upper limits. Upper 
limits were measured using a 10" aperture and a local back- 
ground extracted to match roughly the amount of diffuse and 
scattered emission thought to reside within the source aper- 
ture. Count rates were converted to fluxes with xspec using 
the final ROSAT calibration products and our adopted best-fit 
model. Adopting similar models does not have a significant 
impact on our derived flux upper limits. 

2.3.4. Swift 

The Circinus Galaxy was observed once with the Swift 
X-ray Telescope (XRT) on 2007 March 23 for 7.4 ks. 
Processed data were retrieved from Swift Archive, and 
analysis was performed using FTOOLS and custom soft- 
ware. SN 1996cr is visibly separated from the nucleus 
(25"separation) and diffuse circumnuclear emission and ap- 
pears to be only marginally contaminated. A spectrum of 
SN 1996cr was extracted from the event list using a 10 pixel 
(23 ."5) radius circular aperture, masked to exclude emission 
from the nucleus and diffuse halo to as large an extent as pos- 
sible. Additionally, a background was extracted in a man- 
ner similar to that done for the Chandra observations so as to 
match roughly the amount of diffuse and scattered emission 
believed to be present within the source aperture. Low photon 
statistics, however, limited the fidelity of the background sub- 
traction. We generated calibration products for the spectrum 
and fit it within XSPEC using our best-fit model to estimate 
the flux. 

2.3.5. ASCA 

The Circinus Galaxy was observed once with ASCA using 
the Gas Imaging Spectrometers (GIS) and SoUd State Imag- 
ing Spectrometers (SIS) on 1995-02-14for 61.1 ks. Processed 
data were retrieved from HEASARC, providing w23-36 ks of 
usable data. Analysis was performed using FTOOLS and cus- 
tom software. The target appears unresolved with ASCA and 
was treated as a point source for extraction purposes. Spec- 
tra for the entire Circinus Galaxy were extracted from a 2.'5 
radius circular aperture, while backgrounds were taken from 
blank-sky observations at the same position and with identical 
screening criteria. Contamination from a neighboring source 
5' away was minimal. Total fluxes were measured via simul- 
taneous spectral fitting of the GIS and SIS data and upper 
limits for SN 1996cr were estimated by subtracting the flux 
estimated from the XMM-Newton data as detailed above. 

2.3.6. BeppoSAX 

The Circinus Galaxy was observed on two occasions with 
BeppoSAX using the Low Energy Concentrator Spectrom- 
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TABLE 5 
X-RAY Observations 



Instrument 


Date 


Exp. 


keV 


^2-8 keV 


0.5-2 keV 


2-8 keV 


Comment 


ASCA 


1 nnc AO 1 A 
1995-02-14 


61.0 


< 5' 


< 142' 


< 19^ 


< 152^ 




ROMT HRl 


1995-09-14 


4.1 


< 6 




< 23 






KOSAi HKl 


iyyo-(ji-i6 


1.1 


< 23 




< 93 






ROSA I HKl 


1996-09-10 


1.8 


< 11 




< 44 






ROMT HRl 


1997-03-03 


26.4 


< 2 




< 8 






ROSAT HRl 


1997-08-17 


45.9 


< 1 


— 


< 5 


— 




BeppoSAX 


1995-03-13 


85.2 


< 14t 


< 154"!" 


< 51 ' 


< 165^ 




Chandra ACIS-S 


2000-01-16 


1.0 


13±3 


142 ±35 


62 ±16 


142 ±35 


21% pile-up 


Chandra ACIS-S 


2000-03-14 


4.9 


14±3 


133±15 


51 ±10 


140±16 


2% pile-up 


Chandra ACIS-S 


2000-03-14 


23.1 


10±2 


120 ±20 


45±9 


120 ±20 


20% pile-up 


Chandra HETGS 


2000-06-15 


67.1 


12±2 


108 ±9 


43±9 


118±9 


3% pile-up 


BeppoSAX 


2001-01-07 


27.2 


< 80t 


<614t 


<292t 


< 620t 


Chandra ACIS-S 


2001-05-02 


4.4 


12 ±2 


136 ±20 


53±11 


136 ±20 


21% pile-up 


XMM-Newton 


2001-08-06 


85.5 


16±1 


144 ±5 


61±2 


154 ±5 


Chandra HETGS 


2004-06-02 


55.0 


23 ±3 


186± 11 


88 ±12 


204±11 


3% pile-up 


Chandra HETGS 


2004-11-28 


59.0 


24±3 


196 ±12 


93±11 


215 ±13 


3% pile-up 


Swift XRT 


2007-03-25 


8.0 


35±8 


206 ±39 


124 ±28 


222 ±42 



Note. — Column 1: Satellite and instrument. Column 2: Starting date of observation. Column 3: Exposure time in 
ksec. Columns 4 and 5: Flux in the 0.5-2 keV and 2-8 keV bands respectively, determined from the best-tit model in XSPEC 
in units of lO " erg s ' cm"^. For upper limits, we have adopted the spectrum from the closest Chandra observation. For 
instruments unable to resolve SN 1996cr from the bright nucleus of the Circinus Galaxy, denoted by ^ , upper Umit fluxes 
were estimated by subtracting the total flux of the Circinus Galaxy without SN 1996cr from the total unresolved flux (see 
§2.3 for details and caveats). Columns 6 and 7: Absorption-corrected flux in the 0.5-2 keV and 2-8 keV bands respectively, 
determined from the best-fit model in XSPEC in units of 10"'^ erg s"' cm"^. Column 8: Comments. 



eter (LECS) and Medium Energy Concentrator Spectrome- 
ter (MECS) on 1998-03-03 for 83.9/138.0 ks and on 2001- 
01-07 for 26.9/52.1 ks. Processed data were retrieved from 
HEASARC, providing 83.7/137.2 ks and 26.8/51.7 ks of 
usable data, respectively. Analysis was performed using 
FTOOLS and custom software. The target appears unresolved 
with BeppoSAX and was treated as a point source for extrac- 
tion purposes. LECS and MECS spectra for the entire Circi- 
nus Galaxy were extracted from a radius circular aper- 
tures (limited by the availability of calibration files), while 
backgrounds were taken from blank-sky observations at the 
same position and with identical screening criteria. Contami- 
nation from a neighboring source 5' away was minimal. Total 
fluxes were measured via simultaneous spectral fitting of the 
LECS and MECS data, and upper limits for SN 1996cr were 
estimated by subtracting the flux estimated from the XMM- 
Newton data as detailed above. 

2.4. Archival and Proposed Radio Data 
2.4.1. Australia Telescope Compact Array 

Radio observations of the Circinus Galaxy spanning 1995 
to 2006 were retrieved from the Austraha Telescope Com- 
pact Array (ATCA) archive'^ and are presented in Table 6 and 
Fig. 6. The data were reduced with MIRIAD (v4.0.5) following 
the procedures outHned in the ATNF Miriad User Manual.'** 

The "primary" flux calibrator is 1934-638, which is as- 
sumed to be constant in time with flux densities of 14.94, 
11.60, 5.83, 2.84, and 1.03 Jy at 20, 13, 6, 3, 1 cm, respec- 
tively. As noted in Table 6, there are a few instances in which 
1938-638 is not included as part of the dataset, although the 
observations are initially "primed" using this source; for such 
cases, fluxes of the "secondary" cahbrators were compared 
to their interpolated historical values from calibration mon- 
itoring efforts'' and were fixed to those values when they 

http://atoa.atnf.csiro.au/ 
' httpV/www.atnf.csrro.au/computing/software/mrriad/ 
http://www.narrabri.atnf.csiro.au/calibrators/ 



Strayed by more than ~ 20% from them. The secondary cal- 
ibrators, which served as the gain, bandpass, and phase cal- 
ibrators for the Circinus Galaxy, varied from observation to 
observation and include sources 1921-293 (bandpass, alter- 
nate 1 cm flux), 1329-665 (bandpass), 1549-790 (bandpass), 
1414-59 (phase), and 1236-684 (phase). Both source and 
calibration data were flagged to remove bad time intervals and 
channels, time intervals strongly affected by interference, and 
obvious emission and absorption lines. The calibration tables 
were applied to the Circinus Galaxy data, from which decon- 
volved, primary-beam-corrected images were made. On rare 
occasions, a phase self-calibration was required. Integrated 
flux densities (determined from IMFIT) or 3(t upper Umits (de- 
termined from imstat) were then determined for SN 1996cr. 
These values are provided in Table 6. There appears to be con- 
siderable dispersion in 22.5 GHz flux densities, which were 
mostly acquired as part of larger Australia Long Baseline Ar- 
ray observations. We lack comparison data at other wave- 
lengths to determine if this is simply additional calibration 
uncertainty or potentially real flux variations. However, given 
that there appears to be dispersion even between observations 
taken simultaneously, we caution that the dispersion is not in- 
trinsic, but rather due to systematic errors. 

We note that the only other strong point source in the field 
of view was the nucleus of the Circinus Galaxy, which is 
likely due to emission from the Compton-thick AGN or a 
compact circumnuclear starburst. As a check on the flux of 
SN 1996cr, we extracted flux measurements from the point- 
like nucleus in a manner identical to SN 1996cr. While the 
nucleus could vary intrinsically and may additionally be con- 
taminated at long wavelengths by diffuse emission from the 
extended disk, it still provides us with a secondary estimate 
of any systematic error associated with a given observation. 
We find that the majority of measurements of the nucleus he 
between 10-30% of its mean total flux in each band, although 
a few observations (mainly at 22.5 GHz where contamination 
from water maser emission is possible) vary by a factor of up 
to four; we have adopted larger systematic errors for these ob- 
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TABLE 6 
ATCA RADIO Observations 



Obs. Date 


Obs. ID 


Array 


K-Band 


X-Band 


C-Band 


S-Band 


L-Band 


Comments 












V 




V 


Si, 


V 




V 








C204 


375 


— 


— 


8640 


>0.7 


4800 


>1.3 


— 


— 


— 


— 


lllUsolt, iiClU d. 




C204 


375 


— 


— 


8640 


>0.7 


4800 


>0.9 


— 


— 


— 


— 




1 yyjKu^yj 


C363 


1 5A 


— 


— 


— 


— 


— 


— 


— 


— 


1418 


>6.6 




1 995-04- 1 fi 


C418 


6C 


— 


— 


— 


— 


— 


— 


2768 


>1.1 


1418 


>1.2 






C204 


375 


— 


— 


8640 


>0.7 


4800 


>1.3 


— 


— 


— 


— 






C204 


375 


— 


— 


8640 


>0.7 


4800 


>0.6 


— 


— 


— 


— 


mosaic field "b" 




C204 


375 






8640 


>0.7 


4800 


>1.5 










mr^cQir* fi^^lH 


1995-06-07 


C204 


375 








8640 


>1.5 


4800 


>0.5 














mo^nif fiplfl "Vi" 


1995-07-27 


C204 


6C 














2368 


>1.3 


1376 


>3.4 






CT09 


750B 


















1418 


>17 






C466 


750D 


















1664 


>1.5 




IQQfi-OR-ni 


C363 


375 


— 


— 


8640 


0.5±0.1 


4800 


>1.2 


— 


— 


— 


— 




lQQfi-19-1 ^ 

1 77U 1 Zi- 1 J 


C586 


375 


— 


— 


8512 


2.2±0.3 


4928 


1.8±0.2 


— 


— 


— 


— 


iiujaoic iiciu «. 


lQQfi-19-1 S 

1 77U 1 ^ 1 J 


C586 


375 


— 


— 


8512 


2.0±0.2 


4928 


1.8±0.2 


— 


— 


— 


— 


IIlUsalL ilCiU U 


1 7"U 1 Z,-VJ J 


C505 


375 














2368 


>0.4 


1384 


>0.7 




1QQ7 Of^ 17 


VI 00 


6A 


— 


— 


8425 


14±1 


— 


— 


— 


— 


— 


— 




1 QQ7 OfS 1 7 

vy J I 1 / 


V 1 00 


6A 


— 


— 


8425 


15±3 


— 


— 


— 


— 


— 


— 


apCt-LIcU. IIIIC IIIUUC 


lyy 1 iz. LJj 


V099 


6C 


— 


— 


— 


— 


— 


— 


2268 


22±3 


— 


— 


1QQ7-1 7-0'^ 


V099 


6C 


— 


— 


— 


— 


— 


— 


2268 


37±4 


— 


— 


spectral-line mode 


1997-12-31 


C694 


6C 


— 


— 


8500 


96±15 


4800 


119±12 


— 


— 


— 


— 


O.J vJ n /- CAU.d.pUlalCU 
IIUIII D.W vjnz, Uald. 


1 7 jy Kfy.) 1 7 


C788 


375 


















1384 


89±3 




1 QQQ OfS 1 Q 

1 777 V/U 1 7 


C788 


375 


— 


— 














1418 


85±3 


aJJCt^LIcU. IIIIC IIIUUC 


^yjKjyj 1 7 


V137 


6B 


22190 


76±21 


















9000-0(1-19 


V137 


6B 


22203 


72±20 




















2000-0fi-19 


V137 


6B 


22215 


69±21 






















V137 


1 5G 


22172 


63±10 




















9009-07-1 R 


V137 


1 5G 


22216 


71±10 




















9009-07-1 S 


V137 


1.5G 


22224 


70±10 




















900'^-0Q-07 


C1224 


EW367 


23659 


63d=10 


















nViflftf* fill finlv fliiTT fnffpH 


200^-1 1-0^ 


CI 049 


H214 


16960 


SOitlO 


















■JOO'^-1 1-0"? 


C1049 


H214 


19008 


77±10 




















200^-1 1-04 


C1049 


H214 


21056 


79±10 




















200"^-! 1-04 


C1049 


H214 


22796 


73±10 




















200'?- 1 1-OS 


CI 049 


1.5D 


— 


— 


8256 


163±10 


5056 


227±11 


— 


— 


— 


— 




9004 04 09 


CI 424 


6A 




— 


8640 


162±8 


4800 


240±10 


2368 


364±16 


1384 


481±19 




9004-0S-90 
ZAjy}^ yfj ZA) 




1 5B 


22088 


63±15 


















JJIla&C Col UIlIJ'j ilUA lUlCCU 


2004-06-13 


V176 


750D 


22209 


57±13 


















2004-06-13 


V176 


750D 


22195 


59±11 




















2004-08-01 


V176B 


H168 


22209 


66±6 




















2004-08-01 


V176B 


H168 


22195 


66±6 




















2005-03-14 


CI 368 


H214 


22193 


51±15 




















2005-04-20 


C1368 


750A 


22194 


68±6 


















phase cal only 


2005-11-13 


C1049 


1.5C 






8640 


159±5 


4800 


246±12 










2006-04-29 


C1049 


H214 


18752 


77±7 




















2006-04-29 


C1049 


H214 


21056 


67±6 




















2006-06-21 


C1049 


1.5D 






8640 


154±11 


4800 


261±10 












2006-09-14 


C1341 


H75 


















1418 


505±20 




2007-06-24 


VX013A 




22316 


67±6 





















Note. — Column 1: Starting date of observation. Column 2: Observing program ID. Column 3: Array configuration. Columns 4, 6, 8, 10, and 12: Mean frequency of observed band in 
units of MHz. Columns 5, 7, 9, 11, and 13: Integrated flux densities (determined from IMFIT) or 3a upper limits (determined from IMSTAT) in units of mjy. Errors include both statistical 
and systematic terms. The systematic error is estimated from the ratio of the measured calibrator fluxes over its estimated historical value based on monitored light curves (typical variance 
was 5-20%). Column 14: Comments. Observations taken as mosaics, in spectral-line mode, or with limited calibration sources are duly noted. 



servations as a result. Thus, aside from the 22.5 GHz data, we 
are generally confident in our measured fluxes for SN 1996cr. 

2.4.2. Australia Long Baseline Array 

We carried out Very Long Baseline Interferometry (VLBI) 
observations of SN 1996cr on 2007-06-24 using the follow- 
ing telescopes of the Australian Long BaseUne Array: ATCA, 
Mopra, Parkes, Hobart, and Ceduna.^" The observations were 
phase-referenced to the nearby calibrator PMN J 1355-6326. 
We observed both senses of circular polarization at a fre- 
quency of 22 GHz with a total bandwidth of 65 MHz, and 
for a total time of ~9 h. We report here only our prelimi- 
nary results and note that a fuU description is forthcoming (N. 

See http://www.atirf.csiro.au/vlbi/ for details. 



Bartel et al., in preparation). SN 1996cr was detected only on 
the three baselines involving the Parkes, ATCA and Mopra 
telescopes. The flux of SN 1996cr is consistent with previ- 
ous ATCA measurements, while its extent appears resolved 
by the longest baselines at > 3(7. Although Gaussian and disk 
models provide adequate fits, we adopt a model consisting of 
an optically thin spherical shell, as was found to be appro- 
priate for SN 1993J (see, e.g., Bietenholz et al. 2001; Bartel 
et al. 2002). For a ratio of the outer to the inner shell radius 
of 1.25, we find an outer angular radius of ~5 mas, which 
corresponds to ~ 2.8 x 10'^ cm for an assumed distance of 
3.8 Mpc. Modehng a smaller outer to inner shell ratio Uke 
1.10 could perhaps narrow the angular radius by ^5%, while 
using a filled-center model (equivalent to an inner shell radius 
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of 0) might increase the angular radius by ^20%. 

2.5. Archival Gamma-Ray Burst Data 

Among the numerous GRB-associated SNe that now ex- 
ist, nearly all have been identified as type Ibc's, which sug- 
gests that they arise from Wolf-Rayet progenitors. Although 
SN1996cr is identified as a type Iln SN at late times, it 
is worth searching for any evidence of temporal and spatial 
gamma-ray emission, as its original type may have been quite 
different. Additionally, there are at least two events in which 
a good observational case can be made for the association of 
a GRB with a type Iln SN: GRB 011121/SN2001ke (Gar- 
navich et al. 2003) and GRB 970514/SN 1997cy (Germany 
et al. 2000). We find eight BATSE GRB detections over the 
lifetime of the Compton Gamma-Ray Observatory (Meegan 
et al. 1996) which include the Circinus Galaxy within their 3cr 
error regions. One of these GRBs, 4B 960202, appears to lie 
within the explosion date constraints of SN 1996cr, suggest- 
ing a tentative association. Importantly, however, 4B 960202 
was detected by multiple Gamma-ray observatories and has a 
dramatically smaller error estimate via triangulation from the 
Interplanetary Network (IPN; Hurley et al. 1999, 2005). The 
full IPN 3(7 confidence annulus is only «0.°33 wide and lies 
> 1 .°3 from the SN at its closest tangential point, thus strongly 
ruUng out any identification with SN 1996cr. 

3. ISOLATION OF SN 1996cr'S EXPLOSION DATE 

Tables 3-6 give a detailed summary of the observations rel- 
evant for determining the explosion date of SN 1996cr. While 
SN 1996cr was discovered in the X-ray band, the archival 
X-ray data available are too limited and subject to potential 
geometrical and structural effects (e.g., early absorption, cav- 
ities) to narrow down the explosion date to better than an 
~ 10 yr window between 1990-2000. Although SN1996cr 
is relatively weU-sampled at radio frequencies, those observa- 
tions similarly suffer from the same early-time effects, thus 
limiting their utility to date the SN. Fig. 6 demonstrates that 
the radio emission begins to sharply "turn on" in mid- 1996, 
signaling that SN 1996cr must have exploded sometime prior 
to this. Several narrow-band optical images from the AAT 
fortuitously bracket the SN explosion date between 1995-02- 
28 and 1996-03-16, allowing us to isolate the SN explosion 
to within a year Final confirmation of SN 1996cr comes in 
the form of the optical spectrum, which signals many telltale 
features of a SN embedded in dense circumsteUar material 
(CSM). 

We note that SN 1996cr's position coincides with a small, 
powerful H II region. However, it is clear from a comparison 
of Figs. 2 and 3 that much of the Ha is due to the unresolved 
SN-CSM interaction (~70%; see also Table 2) and thus is 
Ukely to be associated with photoionized wind material sur- 
rounding the progenitor (e.g., Fesen et al. 1999). Contamina- 
tion is likely to be negligible for other wavelengths. 

The tightest direct temporal constraints arise from compar- 
ing the R-b&nd image taken on 1995-02-28 with the narrow- 
band [SII] image taken on 1996-03-16. As such, we cannot 
constrain the spectral form of SN 1996cr at these early times 
(i.e., whether it is line- or continuum-dominated), which make 
conversions between various bands somewhat uncertain. For 
simphcity, we adopt here an intrinsically flat continuum ab- 
sorbed by E{B — V) = 2.2 mags at early times. This is Ukely 
conservative, considering typical early SNe are intrinsically 
quite blue. The extinction was determined using the three 
narrow-band discovery-image magnitude constraints and is 



R:!0.8 mags more than a standard conversion of our best-fitted 
A^H = 8.4 X 10^' cm~^ from the late-time X-ray spectra; a fac- 
tor of ^2 variance is commonly seen between extinction and 
A^H in our own Galaxy (e.g., Burstein & Heiles 1978), so this 
discrepancy is probably not necessarily reason for concern, al- 
though it could be indicative of further early-time absorption. 
We note that this model should provide adequate conversions 
to broad-band magnitudes since early-time SNe spectra are 
often dominated by strong, broad lines that blend to form a 
pseudo-continuum. 

Using SYNPHOT in IRAF to convert from the narrow AAT 
He II, |0 III], and [S II] bands, we find an equivalent discov- 
ery magnitude of V « 17.8. At the distance of the Circinus 
Galaxy, this equates to an uncorrected absolute magnitude of 
My « —10.2, or Mvc ^ —15 when corrected for the minimum 
expected extinction. This provides a strong lower limit since 
the SN likely went off prior to the discovery images and the 
apparent extinction estimated from the discovery magnitudes 
is substantially higher. 

4. CHARACTERIZING SN 1996CR 

The basic observational picture for core-collapse SNe starts 
at shock breakout (e.g.. Chevalier & Fransson 1994, here- 
after CF94), where the blast wave emerges from the progen- 
itor with a typical velocity of > 10"* km s"' and propagates 
into the CSM formed by the pre-SN steUar wind. The in- 
teraction between the ejecta in the forward shock and the 
CSM leads to the formation of a reverse shock, which trav- 
els back into the expanding ejecta — at least relative to the 
outward expanding contact discontinuity that separates the 
shocked ejecta and the shocked CSM. Fig. 7 provides a ba- 
sic one-dimensional picture of the overall scenario. The 
full three-dimensional structure could be much more com- 
plicated. For instance, the contact discontinuity is suscepti- 
ble to Rayleigh-Taylor instabilities, which are only resolved 
in multi-dimensions (e.g., Chevalier et al. 1992; Dwarkadas 
2000). Likewise, asymmetries could exist in the CSM or SN 
explosion itself due to the nature of the progenitor or the ex- 
istence of a binary companion (e.g., SN 1979C, SN 1987A; 
McCray 1993; Crotts & Heathcote 2000; Monies et al. 2000; 
Sugerman et al. 2005). For a typical slow (v„ ~ 10 km s"'), 
dense (M ~ 10^-10"^ Mq yr"') constant progenitor wind, 
the forward shock/CSM interaction produces a very hot shell 
of shocked CSM (T - 10'' K), while the reverse shock/ejecta 
interaction produces a denser, cooler shell of shocked ejecta 
(r 10^ K). The latter interaction has a much higher emis- 
sion measure and generates copious far-UV and X-ray emis- 
sion which subsequently photoionizes a broad inner ejecta re- 
gion from which many of the optical lines are believed to orig- 
inate. When the CSM density is high and/or the power-law 
density distribution of the ejecta (oc R~") is steep, the reverse 
shock will remain radiative and a relatively thin, cool, dense, 
partially absorbing shell (CDS) should form between the re- 
verse shock and the contact discontinuity. UV and X-ray ra- 
diation associated with the original shock breakout (detected 
for flie first time from SN 2008D with Swift; Soderberg et al. 
2008), as well as radiation from the ongoing forward and re- 
verse shocks, may additionally photoionize the outer CSM. 
When the CSM is H-rich, this ionization should give rise 
to the strong, narrow-line H emission which is thought to 
epitomize the type Iln class (e.g., Schlegel 1990). To place 
SN 1996cr within the context of the above scenario, we ex- 
amine our temporal and spectral data in detail. 
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Fig. 6. — The solid circles denote the light curves for V-band (top panel), radio (middle panel), and absorption-con'ected X-ray (bottom panel) data. Fluxes 
are shown on the left, while luminosities are on the right. The light grey shaded region on the left denotes the current explosion window. For the V-band, we 
also show a typical V-band Hght curve based on five type Iln SNe (dashed blue line). There is a large range in properties at early times, for which we show only 
the upper and lower extremes expected; SN 1996cr's maximum light could lie anywhere in between. We have normalized to fit the initial discovery magnitudes 
and adopted a late-time decline rate of AV = 0.0035 per day to fit the HST data. None of the magnitudes has been corrected for extinction, which is likely to 
be in excess of A\j « 4.5. The dark grey region denotes the range of parameter space where our adopted SN light curve is not excluded. The radio data provide 
relatively robust constraints on emission at early times throughout the explosion window, and only peak up in late 1996. The radio emission is fit with a modified 
ejecta-CSM interaction model (dashed lines, see § 4.3), which deviates from the standard model in two important ways. First, the unusual sharp rise and variable 
spectral index cannot be fit by free-free absorption alone, and requires a >100 fold increase in radio luminosity as well. Second, the late-time slope remains flat 
up to the current epoch, rather than exhibiting the characteristic rollover seen in most RSNe. The late-time X-ray data show a relatively strong and significant 
rise between 2000-2007 both in the 0.5-2 keV and 2-8 keV bands. This is very atypical of X-ray detected SNe and suggests that the blast wave of SN 1996cr is 
encountering a dense circumstellar shell. The dashed lines show empirical power-law fits to the data with labeled slopes. The early-time upper limits in all bands 
suggest that either very strong absorption played a role early on or more likely the progenitor of SN 1996cr formed a cavity during the final part of its evolution. 
The thin solid curves, the colors of which mirror those describing the SN 1996cr data, show the corresponding X-ray and radio light curves of SN 1987A for 
comparison (Hasinger et al. 1996; Manchester et al. 2002; Park et al. 2006). Only the luminosity and days after explosion axes apply to the SN 1987A data. The 
SN 1987A data have additionally been scaled up by a factor of 1,000 for illustrative purposes. While many of the specifics are different, the general trend of the 
SN 1996cr X-ray and radio Hght curves resembles those of SN 1987A. 
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Fig. 7. — Idealized one-dimensional cartoon of a SN and its shocks, along 
with the stellar wind established CSM, and in the case of SN 1996cr a more 
distant H 11 region the SN is likely to be suiTounded by. The forward shock 
and shocked CSM are thought to give lise to radio, UV, and hard X-ray emis- 
sion (T ~ 10' K), while the reverse shock and shocked ejecta provide the 
bulk of the soft X-ray emission (7 ~ 10^ K). A CDS immediately behind the 
contact discontinuity can potentially absorb the soft X-ray emission, leading 
to a variety of photoionized emission lines. Not to scale. 



4.1. Optical Light Curve 

To assess the early history of SN 1996cr, we compared the 
V-band light curves for five type Iln, as well as 25 additional 
type 11 SNe of varying types, to the optical constraints de- 
tailed in §2 and plotted in Fig. 6."' The three V >15 upper 
limits from R. Evans here provide some additional insight into 
the nature of SN 1996cr's light curve. We note that the light 
curves of type 11 SNe are relatively uniform and can generally 
be divided into two distinct parts: (1) a sharp rise and vari- 
able decline rate between AV = 0.005-0.035 per day within 
the first ^100 days, followed by (2) a nearly universal decline 
from Mv - -14 with a rate of AV = 0.003-0.01 per day. The 
latter is thought to be linked to the decay of ^''Co (Turatto 
et al. 1990; Patat et al. 1994). 

The above parameters define our average template light 
curve, and we use as our pivot point the constraints from 
1996-03-16 to 1996-03-19, which He 0-380 days from maxi- 
mum. As a consistency check, we note that the slope between 
our 1996 and 1999 data points falls roughly within the ex- 
pected range for type 11 SNe; we adopt AV = 0.0035 per day 
to fit our points. We find that R. Evans' constraint on 1995- 
03-31 excludes earlier explosion maximum dates for all but 
the faintest observed light curves, while the limit on 1995- 
08-13 excludes brighter than average observed light curves 
up to ^50 days prior to this date, and the limit on 1996-03- 
03 fails to exclude any additional time. The dark grey re- 

SNe light curves were acquired from 
http://virtual.sai.msu.Su/~ pavlyuk/snlcurve/. 



gion in the upper panel of Fig. 6 shows the range of parame- 
ter space which is still viable, demonstrating that these upper 
limits exclude SN 1996cr from being an extremely luminous 
type 11 SN over ^--^15% of the explosion window. From this 
we argue that SN 1996cr ultimately had an uncorrected abso- 
lute V-band magnitude somewhere between -10.2 and -15.5. 
Once corrected for >4.5 magnitudes of extinction, SN 1996cr 
falls within the typical range for luminous type 11 SNe. At 
late times, SN 1996cr lies well above an extrapolation of the 
average decay rate, and is in fact ^1.5 magnitudes brighter 
in 2006 than in 1999. This fact demonstrates robustly that 
SN 1996cr must be buoyed either by circumstellar interaction 
or unresolved light echoes (e.g., Patat 2005; Patat et al. 2006). 
We discuss these two possibilities further in §4.2. 

Several SNe (particularly of the Iln sub-type) have been 
shown to be strong near-lR emitters, due to reprocessed emis- 
sion from hot dust (e.g., Gerardy et al. 2002; Pozzo et al. 
2004). Whether the hot dust arises from preexisting dust in 
the circumstellar gas or from newly formed dust in the ejecta 
is still unclear. The HST N1C3 narrow-band constraints we 
have demonstrate that any near-lR emission above the earlier 
ESO near-lR imaging is likely to be minimal. Thus if a sig- 
nificant hot-dust component formed, it must have either been 
short-lived (e.g., destroyed by the initial UV flash from the 
SN or rapidly cooled) or occurred after the N1C3 observations 
in 1998. Mid-IR observations will ultimately constrain the 
dust-content of SN 1996cr better, since the bulk of anticipated 
dust is likely to be at temperatures of ^100 K (the dust tem- 
perature in SN 1987 A, for instance, is wl80 K; Dwek et al. 
2008). While the pre-SN optical imaging does not provide 
useful constraints on the progenitor due to the strong extinc- 
tion toward the Circinus Galaxy, the ESO near-infrared data 
does not suffer as much. The 7-band limits in 1994 June are 
strong enough to exclude stars with Mj < -7.5, effectively 
excluding the bright end of the luminous blue variable (LBV) 
distribution. 

4.2. Optical Spectrum 

As Fig. 1 demonstrates, the optical spectrum of SN 1996cr 
exhibits strong narrow emission lines of H and O that typify 
the type Iln SNe class (e.g., Schlegel 1990) superimposed on 
very broad emission complexes. The strong narrow lines sug- 
gest that circumstellar interaction plays an important role in 
the overall emission almost from the onset and, because of 
the strength of this emission, it can often mask the true nature 
of the photospheric emission. For instance, SN 1998S (Fas- 
sia et al. 2001) was spectrally a type Iln and photometrically 
a type IIL, SN2002ic (Hamuy et al. 2003; Deng et al. 2004) 
and SN 2005gj (Aldering et al. 2006; Prieto et al. 2008) were 
spectral "hybrids" showing type lln-like H lines superposed 
on otherwise type la-like spectra, and SN 2001 em (Chugai & 
Chevalier 2006) spectrally evolved from a type Ic to a type Iln 
on the timescale of a few years. 

In the case of SN 1996cr, where we have limited temporal 
coverage during the early phases, we cannot deduce its early 
photometric or spectral type. We must rely solely on the late- 
time SN emission lines in our VLT spectrum to provide diag- 
nostics on ejecta abundances and shock emission processes, 
as well as potential information on the mass-loss history and 
evolutionary status of the SN progenitor (e.g., CF94, Frans- 
son et al. 2002, Chevalier 2005). As Fig. 8 shows, aside from 
the heavy extinction, the spectrum of SN 1996cr bears many 
striking similarities to the classic type Iln SNe, SN 1979C 
(Fesen et al. 1999), SN 1986J (Leibundgut et al. 1991), and 
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Fig. 8. — A comparison of the optical spectrum of SN 1996cr witli those 
for a few other type Iln supernovae: SN1979C (Leibundgut et al. 1991), 
SN1986J (Fesen et al. 1999), SN1988Z (Aretxaga et al. 1999). Spectra 
have been smoothed over a 3 pixel scale to facilitate easier comparisons. 
SN 1996cr has strong, narrow components similar to SN 1988Z as well as 
broad asymmetric O complexes like SN 1979C and SN 1986J. 



SN 1988Z (Aretxaga et al. 1999), and thus can probably be 
considered to be representative of the class as a whole. 

We used the contributed iraf package specfit (Kriss 
1994) to model the various spectral components of SN 1996cr. 
We began by fitting a power-law continuum, absorbed using 
the Galactic extinction curve of Cardelli et al. (1989) with 
Ry fixed at 3.1.^'^ However, since there is some degeneracy 
between the power-law slope and degree of extinction, we 
fixed E(B-V) to a value of 1.8 determined by the narrow- 
line Balmer decrement^^ and then determined a best-fitted 
power-law slope of 3.4. This single powerlaw, however, un- 
derestimates the curvature of the continuum between 5500- 
7000 A by wlO%, and thus conversely leads to overestimates 
in some line fluxes and widths. To fit the continuum more ac- 
curately, we instead adopted a broken power-law model with 
slopes of 2.65 above and 4.47 below a break at 6400A, re- 
spectively, yielding a statistically reUable model of the con- 
tinuum over the entire spectrum. Note that the narrow Balmer 
Unes are thought to arise from the unshocked progenitor wind 

While there is known to be some variance in the extinction curves from 
galaxy to galaxy (e.g., Calzetti et al. 1994), such differences should be rel- 
atively minimal above 4000 A. Thus for simplicity we model the extinction 
using this particular Galactic extinction model since a significant fraction of 
the extinction comes from our own Galaxy. 

Note that the decrement has a slight dependence on temperature and 
virtually none on density. 



and thus can be considered as the outer layer of emission. 
We therefore consider the extinction derived from the Balmer 
decrement to be a lower limit and caution that some interior 
components could suffer considerably higher extinction. We 
do not, however, see strong evidence of P Cygni absorption in 
any of the emission lines, suggesting that further absorption 
in the vicinity of SN 1996cr is probably minimal. This may 
also imply that the bulk of the CSM is fully ionized. We find 
only marginal evidence for Na D absorption (^ 1.5c7), which 
can be used to approximate the overall extinction (e.g., Tu- 
ratto et al. 2003). The strength of this absorption appears too 
small compared to our other estimates, although the spectral 
resolution, signal-to-noise, and potential for contamination 
from complex line emission (see below) in the VLT spectrum 
severely compromise any meaningful constraint. We contend 
that the lack of broad-band blue and near-UV detections, com- 
pared to the strong blue continua typically seen in compara- 
ble SNe (e.g., SN 1979C, SN 1980K, SN 1988Z; Fesen et al. 
1999; Immler & Kuntz 2005), is fully consistent with our es- 
timated extinction parameters and hence is reasonable. 

We proceeded to fit all of the obvious emission lines with 
narrow Gaussian components. The parameters of each line 
were originally fit separately, but eventually the central wave- 
lengths and FWHMs of the H, He, and heavier element lines 
were tied together as denoted in Table 7 to reduce free param- 
eters and improve error estimates; this was justified by the fact 
that the free-parameter values overlapped with the nominal 
fixed ratio values to within errors. We note that there were, 
however, several large broad residual wings around the Ha 
and O lines. Thus we fit blended narrow and broad Gaussians, 
again linking central wavelengths and FWHMs as denoted in 
Table 7. We describe each element separately below. 

We should caution that the late-time spectrum we observe 
for SN 1996cr could potentially be contaminated by early- 
time SN spectral features due to the presence of unresolved 
Ught echoes from dust in the immediate circumstellar envi- 
ronment or in the intervening interstellar medium (e.g., Patat 
2005; Patat et al. 2006). Notably, the spectra of SNe near 
maximum Ught are the most likely to continue on as light 
echoes, and thus we would expect such contamination to be 
isolated primarily to the continuum of SN 1996cr since early- 
time spectra of SNe routinely exhibit either blue continua or 
psuedo-continua comprised of extremely broad emission lines 
which trace the high velocity ejecta (Filippenko 1997). For 
SN 1996cr, we find only distinct emission fines of low- to- 
moderate velocity (i.e., < 5,000 km s~^) on top of a smooth 
continuum, with the only notable exception being the region 
immediately surrounding He A5876. The equivalent widths of 
these emission lines are also comparable to other well-studied 
SNe such as those in Fig. 8, suggesting at least that SN 1996cr 
is not unusual. We additionally included a few early-time 
light curves of typical SNe in our spectral fitting, and found 
that aside from some degeneracy with our continuum model, 
such spectra are not generally compatible. Thus we argue that 
the majority of our line estimates, with the possible excep- 
tion of Ha, are unUkely to be strongly contaminated. In the 
case of Ha, there are rare exceptions, such as the type Iln 
SN 1994W and SN 1994Y, where early spectra display bright, 
relatively narrow H lines on top of blue continua from the 
outset (e.g., Wang et al. 1996; Chugai et al. 2004). In these 
cases, the narrow Unes are interpreted as arising from the SN 
shock transmitted into a dense, clumpy CSM very close in to 
the SN (~10^^ cm). Our radio and X-ray constraints at early 
times suggest this is an unUkely scenario for SN 1996cr (see 
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TABLE 7 

LlNL iDLN I Il lCAriONS 



# 


Line 


Component 


F 


pc 


Ac 


FWHM 


1 


[0 11] A7319 


n 


1.4± 0.5 


49.6±17.2 


7324.9±1.4 


=29 


2 




b 


13.6± 1.0 


496.0±34.8 


7324.9 (=1*) 


=29 


3 




cl 


6.2± 0.5 


208.0±16.4 


7420.1 (=30*) 


=30 


4 




c2 


19.3± 0.7 


746.0±26.s 


7251.1 (=31*) 


=30 


5 




c3 


16.8± 0.5 


6b1.8±2U.7 


7190.5 (=32*) 


=30 


6 


[Ar III] A7 1 36 


n 


1.8± 0.2 


76.4± 9.7 


7136.8 (f) 


=8 


7 


He I A7065 


n 


3.1± 0.3 


138.3±12.7 


7065.7 (=21*) 


=21 


8 


rc^ TTn \ /'t^ i 

[S lij A67J1 


n 


3.6± 0.3 


210.1±18.9 


6732. 6±0. 8 


741 ±24 


9 


ro TTn \ /'t i ^ 

[S lij A6716 


n 


4.0it 0.1 


239. 0± 7.8 


6717.6 (=8*) 


=8 


10 


TT T "\/'"/''70 

He 1 A6678 


n 


l.lit 0.2 


70.0±13.8 


6675.7 (=21*) 


=21 


1 1 


[N II] A6583 


n 


19.9± 0.7 


1310.0±45.4 


6582.8±0.2 


=8 


12 


Ha 


n 


34.5± 0.8 


2306.8±51.8 


6563.0±0.2 


669±18 


13 




b 


9.5± 1.0 


634.4±65.2 


6563.0 (=12) 


4077±398 


14 


rxT TTn \ /TC H o 

[N 11] A6548 


n 


6.4± 0.0 


432. 3 ± 0.0 


6547.8 (=11*) 


=8 


15 


[0 1] A6300 


n 


5.9± 0.4 


486.7±30.7 


6296.2±0.4 


=8 


16 




b 


5.3± 0.6 


434.9±51.3 


6296.2 (=15) 


=29 


17 




cl 


2.1± 0.3 


159.8±26.0 


6376.8 (=30*) 


=30 


18 




c2 


5.8± 0.4 


504.3±35.5 


6231.6 (=31*) 


=30 


19 




c3 


2.9± 0.3 


263.7±29.2 


6179.6 (=32*) 


=30 


20 


[Fe Vll] A6087 


n 


0.5± 0.2 


43.8±16.3 


6087.0 (f) 


=8 


21 


TT T \ 

He 1 A5876 


n 


1.6± 0.2 


1 OT \ ^ A 1 

187.8±24.1 


5873. 9±0. 8 


913±65 


22 


nvT TTn A c^cc 

[N 11] A5755 


n 


1.2± 0.2 


152.3±22.2 


5755.0 (f) 


=8 


23 


LO IJ A5577 


n 


0.2± 0.1 


35.6±15.4 


5579.0 (=28*) 


=8 


24 




b 


0.3it 0.2 


Ar\ \ ^ c 

49.3±28.5 


5579.0 (=28*) 


=29 


25 




cl 


0.3ib 0.2 


48.5±26.1 


5644.5 (=30*) 


=30 


26 




c2 


0.3± 0.1 


Ar\ 1 ^ /"\ 

4y.3±2(J.8 


5516.0 (=31*) 


=30 


27 




c3 


0.3± 0.2 


48.2ib31.0 


5469.9 (=32''') 


=30 


28 








1 f^-i-Q? '^ 

i OJo.O^^Z,. J 


5007 Q-l-O ^ 


=8 


29 




b 


3.5± 0.4 


1099.1±130.3 


5007.9 (=28) 


2795±174 


30 




cl 


2.4± 0.3 


682.5±71.4 


5066.4±1.1 


2049±19 


31 




c2 


4.5± 0.3 


1543.0± 117.1 


4952.6±0.6 


=30 


32 




c3 


1.9± 0.3 


690.5±94.2 


4910.8±0.6 


=30 


33 


[0 IH] A4959 


n 


1.6 (=27*) 


540.8 (=27*) 


4960.0 (=28*) 


=28 


34 




b 


1.1 (=28*) 


362.7 (=28*) 


4960.0 (=28*) 


=29 


35 




cl 


0.7 (=29*) 


225.2 (=29*) 


5018.9!(=30*) 


=30 


36 




c2 


1.1 (=30*) 


509.2 (=30*) 


4905.2 (=31*) 


=30 


37 




c3 


0.6 (=31*) 


227.9 (=31*) 


4866.8 (=32*) 


=30 


38 


H/3 


n 


1.9 (=11*) 


745.0 (=11*) 


4861.1 (=12*) 


=12 


39 


He n A4685 


n 


0.4±0.1 


221.8±59.6 


4684.5 (=21*) 


=20 


40 


H7 


n 


0.4 (=11*) 


350.2 (=11*) 


4340.0 (=12*) 


=12 



Note. — Column 1: Component number. Column 2: Emission line. Column 3: Component: Narrow (n), Broad 
(b), Complex (c#). Column 4: Observed flux in units of 10"^*^' ergs s~' cm~". indicates the flux of the line is 
tied to another line according to the relative ratio of their respective atomic line intensities. Column 5: Extinction- 
corrected flux in units of 10~'^ ergs s~' cm~^. indicates the flux of the line is tied to another line according to 
the relative ratio of their respective atomic line intensities. Column 6: Central wavelength in A. indicates the 
central wavelength of the line is tied to another line according to the ratio of their respective atomic rest wavelengths. 
Column 7: Gaussian Full- Width Half-Maximum in km s~\ 



TABLES 
Line Ratios 





Line Ratio 


"n" 


"b" 


"cl" 


"c2" 


"c3" 


[OIII] 


7(4959 + 5007)//(4363) 


> 5.40 


>3.62 


>5.08 


>2.25 


> 2.28 


[Oil] 


7(3726 + 3729)//(73 19 + 733 1) 


< 703 


<3.02 


< 1.78 


<7.46 


< 1.90 


1/[0 II] 


/(7774)//(7319 + 7331) 


<0.50 


<0.16 


<0.12 


<0.48 


<0.34 


[01] 


/(6300 + 6364)//(5577) 


>3.95 


> 1.90 


>3.14 


> 1.07 


> 1.63 


[SII] 


/(6731)//(6716) 












[SII] 


7(67 1 6 + 673 1 ) //(4068 + 4076) 


>0.31 










[Nil] 


7(6548 + 6583)//(5755) 


J 1 44+9.54 










[Ar III] 


7(7136+7751)/7(5192) 


> 0.84 










He I 
He I 


7(7065)//(5876) 
/(7065)//(6678) 


074+0.78 

1.98+3:?^ 

37+^-^^ 










He I 


/(6678)//(5876) 











Note. — 



Column 1: Line-intensity ratio. Columns 2-6: Fitted component values. 
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Fig. 9. — A section of the spectrum in Fig. 1 showing in detail the fitted 
components in the vicinity of the Ha fine. Obvious emission lines have been 
identified and an empirical model has been fit to the data as described in 
§ 4.2; again the green line shows the overall fit to the spectrum while the red 
(blue) Unes show the individual narrow (broad) components. Although the 
narrow Ha strongly overlaps with the N II doublet, our spectral resolution is 
sufficient to deblend it. 

§4.3 and §4.4). However, since we lack early optical spec- 
tra of SN 1996cr, interpretation of the Ha emission should be 
viewed with some discretion. 

4.2.1. Hydrogen 

We find SN 1996cr exhibits both narrow (669±18 km s"') 
and broad (4077±398 km s"') components of Ha, 
with fluxes of (2.31±0.05)x 10"'^ ergs s"' cm'^ and 
(6.34±0.65)x 10"''* ergs s"' cm"^, respectively. This portion 
of the spectrum is shown in detail in Fig. 9. The narrow line 
appears to be marginally resolved above our estimated instru- 
mental resolution, although this fact should be regarded with 
some caution since the line is strongly blended with N II; we 
await confirmation via high-resolution spectroscopic follow- 
up. 

The flux found in our VLT spectrum is consistent to 
within errors with that measured from the HST F656N fil- 
ter seven years earlier, indicating that the overlapping H II re- 
gion still only contributes at most 30% to the total narrow- 
line flux above and 10% to the broad-line flux. Account- 
ing for this contamination, we find corresponding extinction- 
corrected luminosities of (2.79±0.06)x 10""^ ergs s"' and 
(9.86±0.10)x 10" ergs s"' for the narrow and broad compo- 
nents, respectively. The energy generated in either the nar- 
row or broad Ha lines is significantly above that predicted by 
the decay of ^^Co-^^Fe at this late stage (e.g., Chugai 1991), 
and therefore must instead be related either to mechanical en- 
ergy associated with the SN ejecta-wind interaction or poten- 
tial light echoes. 

The line width and systemic velocity of the narrow com- 
ponent argue for an origin beyond the blast wave, associated 
with either the shocked or unperturbed progenitor wind (e.g., 
Fransson et al. 2002). If the line does arise from a dense, pre- 
SN, H-rich CSM, it must have been photoionized and heated 
by UV and X-ray radiation from the massive progenitor and 
the initial shock breakout of the SN. The fact that the Ha lu- 
minosity has remained relatively constant between the epochs 
of the VLT and HST observations suggests that this region is 
now being sustained by ionizing flux from the current ejecta- 
wind shock zone (see §4.4), although again contaminating 
echoes from prior epochs is a concern here. For compari- 
son, the Ha flux of other type lln SNe such as SN 1988Z, 



SN 1995G, and SN 1995N decreased by at least an order of 
magnitude in ~4 yrs (Aretxaga et al. 1999; Pastorello et al. 
2002; Fransson et al. 2002). The Ha luminosity thus provides 
an important diagnostic on the total ionizing radiation that 
complements the X-ray observations, as this line arises as a 
result of recombination and collisional excitation. If this Ha 
component is associated with CSM interaction and is pow- 
ered by photoionization, then the X-ray luminosity from the 
forward and reverse shocks should roughly be proportional to 
the Ha luminosity following Eq. 3.7 in CF94 such that 

L(Ha) «4.4 X 10"(r/10^K)-'(Lx/10'**') ergs s"'. (1) 

This implies Lx = 7.0 x 10'***(r/10^/r) ergs s"', however, 
which is two orders of magnitude larger than the product of 
observed absorption-corrected X-ray luminosity and temper- 
ature. This would argue for the soft X-rays from the reverse 
shock, which should dominate bolometrically, to perhaps be 
radiative and hence reprocessed by either unshocked low- 
density ejecta or high-density material in the CDS (this possi- 
bility is discussed further in §4.2.2). Alternatively, the narrow 
line could arise from shocks being driven into dense, slow- 
moving clumps associated with the progenitor wind. This 
would more naturally explain the discrepancy in observed and 
estimated ionizing flux, in which case the line is instead pow- 
ered by collisional excitation. 

If the narrow Ha stems from the stellar wind, then the Ha 
luminosity can also be related to the mass of the unshocked 
wind as 

L(Ha) = (1 -5) X 10^''(M_4/v.,i)' '"Is ergs s"' , (2) 

where rig is the radius in units of 10'^ cm, v„, is the wind ve- 
locity in units of 10 km s"', and M_4 is the mass-loss rate in 
units of lO"'' Mq yr-i (Deng et al. 2004; Wang et al. 2004). 
We assume an emissivity based on Case B recombination at 
a temperature of 10,000 K, which is in rough agreement with 
our extinction-corrected Ha/H/3/H7 ratios. Our Ha luminos- 
ity implies a mass-loss rate of ~10~^ v^.j r^^ Mq yr~'. For a 
typical red supergiant lifetime of > 10^ yr, this implies that at 
least a few Mq of hydrogen-rich matter was cast off into the 
CSM. 

The high-velocity component of Ha, on the other hand, 
must be related directly to the blast wave in order to have 
achieved its present dispersion. This component could be as- 
sociated with the reverse shock or CDS behind the reverse 
shock, perhaps from swept-up wind material or hydrogen that 
remained on the surface of the progenitor The broad-line Ha 
velocity is thought to provide a lower limit to the true shock 
velocity. However, an observed disperion of ^4,000 km s"' is 
much lower than the 15,000-40,000 km s"' initially expected 
for a typical core-collapse SN, and implies that the shock in 
SN 1996cr has akeady slowed down substantially due to cir- 
cumstellar interaction. This observed velocity decrease of at 
least >4-10 implies a circumstellar density increase of >16- 
100. This line is almost certainly powered by collisional ex- 
citation associated with shocked ejecta. 

We additionally see significant narrow-line emission in H/3 
and H7, as well as marginally in HS. Broad emission com- 
ponents could also exist for these lines; however, for typical 
line ratios with respect to Ha they would be marginal at best. 
Additionally, such emission would be impossible to deblend 
from the [O III] emission complexes. 
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Fig. 10. — Velocity structure of O emission lines in units of km . The 
VLT spectra in the vicinity of five labeled complexes (black) and their re- 
spective model fits (green showing the overall full component fit, while con- 
tributions from a distinct velocity component are shown in different colors) 
are shown. Each spectrum is offset by a factor of 10 for clarity. The solid 
vertical line denotes the approximate position of the stationary narrow (cyan) 
and broad (magenta) components, while the dashed vertical lines denote the 
expansion velocities of components cl (purple), c2 (red), and c3 (pink); a po- 
tential symmetric component c4 is not evident in the data. Although several 
of the lines are strongly blended, each of the components remains relatively 
distinct. The line components for A5007 and A4959 overlap such that in a 
few instances no single component obviously dominates. 

4.2.2. Oxygen 

We find evidence for one narrow-line and four 
broad-line velocity components in the profiles of the 
[O III]AA5007,4959, [O I]A6300, and [O II]A7319 Unes 
as shown in Fig. 10, with velocities of 3577 ±65 ("cl"), 
66 ± 1 8 ("n", "b"), -3247 ± 36 ("c2"), and -5755 ± 37 km s"' 
("c3"). We note there is also marginal evidence for similar 
multi-component emission complexes associated with the 
[O I]A5577 and [O III]A4363 lines. This is somewhat in 
contrast to the late-time spectra of many type 11 SNe such 
as SN 1993 J, which show evidence for more irregular, 
"clumpy", newly synthesized O complexes (Matheson et al. 
2000). 

A few type 11 SNe have shown evidence for three-peaked 
profiles, which have either been ascribed to inhomogeneous, 
partially absorbed disk- or torus-shaped CSMs (e.g., Gerardy 
et al. 2000; Leonard et al. 2000), or from line formation in 
the CDS (e.g., Fransson et al. 2005). In such cases, how- 
ever, the profiles arise not from the O lines, but from Ha, 
Mg II, or Lya; this may indicate that SN 1996cr is O-rich. 
Nonetheless, the first scenario typically yields two relatively 



symmetric positive and negative velocity peaks indicative of 
emission from an equatorial ring, while a third central peak 
could emerge due to a population of dense clouds embedded 
in a much less dense CSM. In the second scenario, the CDS 
can only lead to such peaks under a very specific set of con- 
ditions such that it is very optically thick to such line emis- 
sion, it is thin enough that the velocity gradient over the shell 
is small compared to the thermal velocity, and finally it has 
a macroscopic velocity larger than its thermal velocity. The 
CDS could be clumpy, lowering the filling factor and allow- 
ing some emission to escape. A final possibility for this ve- 
locity structure could arise from the SN blast wave impacting 
a dense shell produced by a wind-blown bubble. In this case, 
several secondary shocks could separate the CSM into multi- 
ple concentric shells (e.g., Dwarkadas 2005), such that each 
peak may represent a different velocity shell. The inner (red- 
shifted) shells would be preferentially absorbed compared to 
the outer (blue-shifted) shells, just as we observe. 

For SN 1996cr, we see a broad central peak ("b", 
FWHM=2795±174 km s"') and two roughly symmetric pos- 
itive and negative peaks ("cl" and "c2", FWHM= 2049 ± 
19 km s"'). The red peak is always fainter than the blue one, 
implying significant internal extinction. Interestingly, we also 
see a fourth, faster negative component ("c3") which does not 
have an obvious positive counterpart; with our current sen- 
sitivity such a symmetric component (i.e., "c4?") must be 
at least a factor of 2.5 fainter (3(t), implying significant ad- 
ditional internal absorption above that of the lower velocity 
peaks, if it exists. This "c3" component could be indicative 
of further asymmetry in the structure of SN 1996cr. We need 
only look to the nearby SN 1987 A, with its triple ring struc- 
ture, for an example of such potential complexity (McCray 
1993, 2007). While the properties of the CSM of SN 1996cr 
appear to differ dramatically from those of SN 1987 A, the 
spectral features hint at an environment that is equally dy- 
namic and exotic. Another possibility is that the "c3" com- 
ponent represents a light echo from an earlier epoch, where 
the O lines were shifted by a larger velocity and more heavily 
absorbed. We reject this, however since it would require the 
"c3" oxygen features to peak strongly in the blue part of the 
spectrum, while we detect them clearly throughout the spec- 
trum. 

Notably, when we apply our five-component model to the 
O lines, we find that our constraints on [O I] A6363 are consis- 
tent with zero flux if we fit [O I]A6363 and [O I]A6300 sep- 
arately. Given the degeneracy between the potential broad 
components of these two lines though, this only provides 
weak /(A6363)//(A6300) constraints for the various compo- 
nents (n:< 1.0, b:< 0.8, cl< 0.6, c2< 1.5, c3< 0.9). Such 
ratios are all fully consistent with the normal transition prob- 
ability of 0.3. For the following, we fit only A6300 and use 
its flux to estimate the total intensity /(A6300+ A6363) for 
the various components. As such we may incur some small 
additional errors associated with flux contamination between 
different velocity components of the two lines. We also cau- 
tion that the complex which we attribute solely to [O II]A7331 
is likely in fact to be a blend of [O II]AA7319,7331 and 
[Ca II]AA7291,7325. Using our five-component model, how- 
ever, we find that additional lines do not significantly improve 
X^- The relative symmetry and constraints on line centers 
of the individual components imply that other lines do not 
contribute strongly to the overall emission, although notable 
small systematic residuals suggest minor blending is likely. 
The overlapping broad-line O profiles seen in Fig. 10 appear 
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to be overestimating the valleys in between distinct compo- 
nent peaks, indicating that either the profiles are not repre- 
sented well by Gaussians or that P Cygni-Uke absorption may 
be present but masked by the overall complexity of the emis- 
sion. Higher-resolution spectroscopy is needed to place better 
constraints on these potential emission and absorption contri- 
butions. 

The fact that the various components show up in all of the 
O lines strongly argues for a scenario in which they arise 
from physically associated regions within the SN. Since we 
see several ionization states though, it is likely that there is 
strong stratification within these distinct spatial regions (very 
likely due to an extended partially-ionized zone; e.g., CF94) 
and it is only appropriate to employ nebular-to-auroral line ra- 
tios for a given ionic species to estimate physical conditions. 
However, our constraints on AA3726, 3729, A4363, and A5577 
are poor and do not provide physically interesting temperature 
and density limits. Table 8 lists these oxygen ratio constraints. 
Notably, the relative observed strengths of the O lines are 
broadly consistent with the model predictions of CF94 (see 
Figs. 4 and 5), whereby the modeled ejecta is photoionized 
by X-rays from the circumstellar interaction. 

Finally, there appears to be no sign of the O I A7774 re- 
combination Une either at its systemic value or other veloci- 
ties seen in the [O I] lines (although telluric absorption does 
hinder assessment of the blue components). This suggests 
that the narrow O resides in a physically distinct region sep- 
arate from the narrow Ha, as one typically sees a 10"^ ratio. 
The broad, high-velocity O emission is not mirrored in any 
other elements, although there is possibly some complexity to 
the He emission (see below), further arguing that it probably 
arises from an O-rich shell of processed gas associated with 
freely expanding ejecta. 

4.2.3. Helium 

SN 1996cr shows evidence for four prominent, narrow He 
lines (906±101 km s"'): He IA7065, Hel A6678, Hel A5876, 
and He II A4686. The narrow-line width is marginally broader 
than that of Ha, hinting that it may be spatially distinct. The 
He I Une ratios do not strongly constrain the electron tem- 
perature and density of the gas. The emissivities of Porter 
et al. (2005) constrain the He gas to be at T > 10'' K and 
n^, > lO'* cm-\ 

We see substantial broad residual emission around A7065 
and A5876 which currently defies identification. Unhke the 
well-formed velocity components associated with the vari- 
ous O lines, there appears to be no spectral consistency be- 
tween different He lines. For instance, the emission around 
A7065 appears box-like and quite strong relative to the sys- 
temic narrow line, while the emission around A5876 appears 
more spread out and much fainter relative to the systemic nar- 
row line. Thus if this broad emission is associated with He, 
it is likely present in several regions of the SN, and may be 
far more complex than H or O are. There may be additional 
marginal broad features around A4686, while any potential 
features around A6678 are masked by the strong, broad Ha 
emission. 

The narrow He emission likely arises from the same re- 
gions of the progenitor wind as the narrow hydrogen emis- 
sion, whereas the broad unidentified emission, if it comes 
from heUum, could be swept up material or expelled clumps. 
Either way, its complex profile suggests it is much more irreg- 
ularly distributed compared to H. 



4.2.4. Nitrogen 

We find evidence for strong, narrow N emission in the 
spectrum of SN 1996cr in the form of both the nebular 
[N II] AA6583,6548 (which appear as a resolvable blend with 
Ha) and the auroral [N II]A5755 lines. The lines have a 
best-fitted FWHM of 741 ±24 km s~\ The hne-intensity ra- 
tio /(A6583 + A6548)//(A5755) implies an electron density of 

> 1 X 10^ X (104/r)°-5 cm-3 for T < 20,000 K (Osterbrock 
1989; Keenan et al. 2001). The high density suggests that 
clumping of N is likely. As argued for the narrow H emission, 
the majority of the emission must arise from the SN-CSM in- 
teraction given the pre-SN imaging constraints. The implied 
velocity of the lines suggest it is likely to have a circumstel- 
lar origin, while the strength of the emission hints at poten- 
tial nitrogen enrichment, as has been notably seen in several 
other type II SNe such as SN 1979C, SN 1987A, SN 1993J, 
SN 1995N, and SN 1998S (Fransson et al. 1984, 2002, 2005). 

4.2.5. Sulfur 

The spectrum shows strong, narrow [S II]AA6731,6716 
emission lines with a best-fitted FWHM of 741 ±24 km s"\ 
which again clearly must come from the SN-CSM iteraction 
given the pre-SN imaging constraints. This feature is rela- 
tively uncommon in late-time type 11 SNe spectra, although 
it is seen in several comparatively-aged type lln SNe spectra 
such as SN 1988Z and SN 1979C (see Fig. 8). The [S II] hne- 
intensity ratio /(A6731)//(A6716) implies an electron density 
of n, = 400!^™ X (lO'^/T)"-^ cm'^ (Osterbrock 1989). We 
additionally see a hint of the [S II]AA4076,4069 doublet, al- 
though the significance of this feature is marginal and strongly 
affected by our estimated extinction value. Given the strong 
density contrasts, this component should reside in a different 
region from [N II]. 

4.2.6. Argon 

Narrow Ar emission is present in the spec- 
trum of SN1996cr in the form of both the nebular 
[Ar III]AA7751,7135 and the auroral [Ar III] A5 192 
lines. Our lower limit to the line intensity ratio 
/(A5192)/7(A7751 + A7135), however, is too loose to 
place any useful constraint on the electron temperature or 
density of this component (e.g., Keenan et al. 1988). 

4.2.7. Iron 

We find significant high-ionization [Fe VII]AA6087,5721 
lines, but only a marginal [Fe VII]A5158 line. Additionally, 
[Fe X]A6375 may be present, although it is currently lost 
amid the broad [O I] emission. The lack of a strong A5158 
line constraint unfortunately leads to physically uninterest- 
ing electron-density limits (e.g., Keenan & Norrington 1987, 
1991). 

4.2.8. Unidentified 

We also see tentative features at A5840, A5980, A701 1, and 
A7785 which remain unidentified. 

4.3. Radio Light Curve 

An expectation from the interaction between the blast wave 
and the CSM is the production of copious non-thermal syn- 
chrotron emission from relativistic electrons and enhanced 
magnetic fields within the thin shell trailing the forward 
shock. The radio emission should be proportional to the in- 
jection spectrum of electrons and various loss mechanisms, 
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which are effectively determined by the composition and den- 
sity of the CSM. At early times, however, the radio emis- 
sion can be absorbed via free-free (FF) or synchrotron self- 
absorption (SSA) depending on the magnetic field strength 
and CSM density and structure. Then as the shock wave 
overtakes more and more of the CSM, the radio emission be- 
comes progressively less absorbed, leading to a characteristic 
frequency-dependent "turn on" first at higher frequencies and 
later at lower ones. 

Following the prescription of Chevalier (1982b) the opti- 
cally thin radio luminosity of this shell can be described by 



U oc AttR^AR K B 



(7+l)/2 j^(7-l)/2 



(3) 



where Ly is the radio luminosity at frequency v, AR is the 
thickness of the synchrotron emitting region at radius R, B 
is the strength of the magnetic field, the distribution of ac- 
celerated particles is assumed to take a power-law form of 
A^(£') = KE~^, and t is the absorption opacity. Both B and K 
scale with the thermal pressure P, which itself scales with the 
CSM density pcsM of the synchrotron emitting region. The 
opacity t likewise scales with pcsM, although in the case of 
FF absorption it is proportional to the integrated density along 
the line of sight. We assume here that all absorbing media 
are purely thermal, ionized hydrogen with opacities propor- 
tional to We contend that the FF opacity dominates 
over the SSA opacity given that (1) the shock velocity re- 
quired by SSA, as implied by the radio luminosity and time of 
peak emission (see, e.g.. Chevalier et al. 2006, and references 
therein), is many times lower than either the optically inferred 
late-time shock velocity or the expected early-time shock ve- 
locity, and (2) the absorption does not deviate substantially 
from T oc ' where it is possible to constrain it, as expected 
for FF. By comparison, SSA should scale as tssa oc v"'^-^, 
which is too steep to fit any of our well-sampled epochs. In 
this case, the FF opacity is given by 



-1.35 



Tff «3.28 X 10" 



-2.1 



EM 
pccm" 



where the emission measure, EM, is given as 



(4) 



poo poo / \ 2 

EM= / N^dr= / (^) dr. (5) 

Jr Jr V AiH'«H / 

Thus the unabsorbed radio luminosity and FF opacity Tff 
afford us two completely independent opportunities to trace 
the CSM density profile. 

With the above in mind, we now examine the radio data for 
SN 1996cr. Fig. 6 demonstrates that there are strong upper 
limits at four observable frequencies from the earliest time of 
explosion up through the end of 1996. These non-detections 
imply either strong early absorption or the presence of a low- 
density cavity. Over the next year, we see a dramatic rise 
in the radio emission, first at higher frequencies and later at 
lower ones. The inverted rise, however, is atypical and cannot 
solely be the result of dwindling absorption, as is common 
in other RSNe. To make this point clearer, we present the 
radio spectrum of SN 1996cr for several epochs in Fig. 11. 
Much of the radio data were taken separately, and thus a 
few of our adopted late-time epochs actually span a several- 
month window; such an approach is validated by the fact that 
strong, rapid radio variabihty is not expected or observed at 
late times. 
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Fig. 11. — Radio spectra of SN1996cr over several epochs. Solid cir- 
cles denote data, while solid and dashed lines denote intrinsic and absorbed 
radio spectra from our best-fitting model. We find dramatic changes in the 
low-frequency slope indicative of time-varying absorption, but relatively little 
change to the high-frequency slope additionally suggesting a strong increase 
in the CSM density. 



From Fig. 11, we see that the data taken during the 2006 
June and 2004 April epochs are only minimally absorbed 
and thus provide a solid constraint both on the intrinsic syn- 
chrotron spectral index (a = -0.79 ± 0.02) and the FF ab- 
sorption beyond the synchrotron-emitting region. Using our 
model radio spectra to guide the eye (described below), we see 
the presence of increasing absorption at the lowest frequen- 
cies as we progress back through the 1999, 1997, and 1996 
epochs. This is a natural consequence of the radio-emitting 
region overtaking the CSM which absorbs it. Importantly, the 
observed spectral index between 4.8 GHz and 8.5 GHz does 
not change dramatically between 1996 and 2006. If we as- 
sume that the intrinsic spectral index is constant, this implies 
that the spectrum here is optically thin. The 8.5 GHz emis- 
sion is thus only modestly affected by FF absorption (^ 20%) 
even at early epochs and can be considered a relatively robust 
tracer of the intrinsic radio luminosity at all times. As such, 
the 8.5 GHz flux density appears to jump by a factor of ~150 
between 1996 August and 1997 December (515 days) and an 
additional factor of rs1.33 between 1997 December and 2004 
April (2384 days). Given that the radio-emitting region is 
spread over relatively large scales, the radio emission should 
vary smoothly in time. Empirically, we find that a broken 
power law, with indices ai = 8.76 ±0.24 and 02 = 0.71 ±0.09 
between the three highly constrained epochs above, provides 
an adequate fit to the data from other epochs aside from 1996 
August (xl = 2.65), although there is clearly some degener- 
acy with other parameters such as a and Tff (see Fig. 1 1). This 
model should at least provide a qualitative understanding of 
the light curve during this overall period. 

The early upper limits and dramatic increase in the radio lu- 
minosity imply a sharp rise in pcsM as well, highlighting the 
possible transition from a fast, sparse stellar wind to a slow, 
dense one. Such transitions typically lead to the formation of 
a wind-blown bubble (e.g.. Weaver et al. 1977; Garcia-Segura 
et al. 1996a,b), and numerous researchers have explored the 
subsequent interaction between the complex CSM associated 
with a bubble and the SN blast wave (e.g.. Chevalier & Liang 
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1989; Tenorio-Tagle et al. 1990, 1991; Dwarkadas 2005). 
Quantitative constraints on the density are difficult to obtain, 
however, since the radio luminosity cannot be related to the 
CSM density in the standard manner (i.e., self-similar solu- 
tions such thatLy oc Pcsm~Pcsm for ^=0.79; Chevalier 1996) 
and requires hydrodynamical simulations to account for the 
interaction properly. 

Our constraints on ts, on the other hand, are more straight- 
forward to calculate. Using the Levenberg-Marquardt least- 
squares method we directly constrain tq = rsiv/l GHz)^ ' to 
be 29.3±6.1, 15.8±1.1, 4.1±0.6, 0.69±0.10, and 0.57±0.11 
for the 1996 December, 1997 December, 1999 June, 2004 
April, and 2006 June epochs, respectively. Unfortunately, the 
1996 August and 1997 June epochs only have 8.5 GHz mea- 
surements and do not provide strong constraints, although we 
are encouraged by the fact that the latter is consistent with our 
adopted light-curve model. For the five epochs where we have 
good spectral constraints, we solve for pcsM,sheii explicitly in 
four contiguous shells such that 

io^k; v/; mdrJ 

(6) 

where we have assumed /xh = 0.6 and Te =10,000 K. The 
inner and outer radii r" and in pc, are taken from 
Vjf ;^'* we note that the outermost radius is constrained to be 
« 2.8 X 10'^ cm on 2007-06-24 from our VLBI observation. 
The to-be-determined r-dependence of pcsM,r is encompassed 
by f(r). The form of /(r) is iteratively determined by ensuring 
that it is consistent with the slope obtained between adjacent 

shells of PcSM.shell- 

We next turn to the evolution of the shock velocity, Vj, 
which we need in order to determine the shell radii. From 
our optical spectrum, the ejecta (as traced by the O lines) 
and broad Ha line appear to be expanding at a variety of 
velocities between w2000-5800 km s"', with the symmetric 
red and blueshifted velocities of ?s3250-3550 km s"' stand- 
ing out as the most likely values tracing the blast wave. We 
thus adopt 3400 km s"' for the current shock speed. Clearly 
after the shock impacts a high-density region, its velocity will 
be reduced considerably. During its subsequent evolution in 
such a region, however, its velocity should evolve strongly, 
as can be seen from radio observations of SN 1987A (Manch- 
ester et al. 2002; Gaensler et al. 2007) and simulations thereof 
(Dwarkadas 2007b). Therefore, without any further informa- 
tion, a first approximation is to assume a constant velocity 
comparable to that derived from the optical spectrum. Inter- 
estingly, the high X-ray temperature, if it is in thermal equi- 
librium, also yields an equivalent velocity of f«3400 km s~\ 
although here the velocity is almost certainly related to the 
reverse rather than forward shock and thus may simply be co- 
incidence. 

We therefore adopt an evolution for the blast wave velocity 
as follows. The shock began with a high initial speed which 
it maintained for 1-2 yrs. By 1996 August it encountered a 
dense shell of wind-swept material, and the shock velocity 

dropped rapidly by a factor proportional to the square root 
of the density gradient it encountered. This density gradi- 

While a self-similar solution is unlikely to hold under our current phys- 
ical conditions, it is instructive to note that in such a scenario, the radii and 
velocity would evolve as r = (v,/m)< and v., oc r"\ where m will range be- 
tween 0.85-0.95 for most cases. This would lead to estimated shell radii 
only ~ 20% larger and should highlight the relative quality of our radial con- 
straints. 



ent also initiated a rise in the radio luminosity by a factor of 
> 200. Beyond 1996 August, we assume the velocity quickly 
arrived at a constant value of ~3400 km s"' . These assump- 
tions allow us to determine approximate radii of 0.054 pc, 
0.055 pc, 0.059 pc, 0.064 pc, 0.081 pc, and 0.089 pc for the 
1996 August, 1996 December, 1997 December, 1999 June, 
2004 April, and 2006 June epochs, respectively. 

If we take the radius on 1996 August as our shell impact 
date, then, depending on the 380 day ambiguity in the explo- 
sion date, the average initial shock velocity prior to this is con- 
strained to be ~37,000-l 39,000 km s^^ (0.12c-0.46c). The 
lower value is in line with theoretical and observational expec- 
tations (e.g.. Chevalier et al. 2006), while the upper value is 
probably unrealistic. Such a high average shock speed would 
be difficult to maintain for more than several days in a typical 
or even underdense CSM environment. If the explosion did 
come later, then one or more of our assumptions above may 
be incorrect. The late-time velocity of SN 1996cr could be 
larger than our adopted value: the largest observed velocity 
in the optical spectrum is f«5,800 km s~\ for instance, and 
provides a plausible upper limit to the current shock speed; 
adopting this value yields initial velocity constraints a factor 
of ~3 lower. Alternatively, the transition between initial and 
final velocities could have been gradual, taking many weeks 
to years to fully drop from one extreme to the other: the 
proposed CSM shell of SN 1996cr could have asymmetries 
akin to those seen in SN 1987A's ring, for instance, which are 
hkely responsible for the onset of hotspots (e.g., Sugerman 
et al. 2002); a long-lived velocity transition could again lower 
initial velocity constraints by a factor of a few. Finally, the 
VLBI radius constraint could be too large, due either to our 
adopted ring model or our adopted distance to SN 1996cr: we 
note that uncertainties in our adopted VLBI model are esti- 
mated to be only ~5-20%, while our adopted distance could 
vary by a factor of <2.7 at most. Despite these potential 
problems, though, it is encouraging that our overall velocity 
model, assuming the lower initial shock velocity at least, is re- 
markably similar to the actual time evolution of the shock ve- 
locity in SN 1987 A as deduced from radio observations (e.g., 
Manchester et al. 2002). In SN 1987A, the shock expanded 
at ~35,000 km s~^ for its initial three years, before impacting 
with the H II region and slowing to '--^4500 km s"' . 

Using Equation 6, we can now place constraints on pcsM.r 
in each shell of material overtaken between the five epochs 
as shown in Fig. 12. Since we do not know how far out the 
remainder of the free-free absorbing material extends into the 
CSM, we do not include the last tq value in this analysis. Ha 
observations from the Southern H-Alpha Sky Survey Atlas 
(Gaustad et al. 200 1 ) provide constraints on the emission mea- 
sure toward the Circinus Galaxy of « 10,000-25,000 pc cm"^, 
or equivalently tq = (0.0033-0.0082)(7;/10'' K)''^-^^. Thus for 
any reasonable temperature of the warm interstellar medium 
(ISM), it appears that much of the remaining free-free absorp- 
tion from our last measurement probably resides close to the 
SN in the CSM itself. Future low-frequency radio measure- 
ments could constrain this potential decline in the absorption. 
Oddly, we find that the /OcsM,r oc profile derived from the 
absorption, where s « 2.5-2.9, seems at odds with the flat 
late-time evolution of the radio luminosity itself, which would 
naively imply i w 1.^^ The strong disconnect between the two 

Radiation losses or strong gradients in the magnetic field could change s 
substantially (e.g., Fransson & Bjomsson 1998). Note, however, than any ad- 
ditional loss mechanism would have to leave the late-time spectrum relatively 
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Fig. 12.— Estimated CSM density profile for SN 1996cr; see §4.3 for full 
details. Radii are determined assuming a constant velocity of 3400 km s' 
since 1996 December, anchored by our VLBl size constraint shown as a 
thick grey line. The CSM densities are estimated from the decreasing free- 
free absorption coefficient assuming a fully-ionized CSM (see 4.3 for de- 
tails), and hence only trace the ionized portion of the CSM at these radii. 
The green curve denotes the density profile obtained assuming pcsM °c f^, 
which is clearly inconsistent with this assumption. The blue points and curve 
show the density profile when modeled in a self-consistent manner, such that 
PCSM <x where i is the iteratively determined slope between subsequent 
data points to be —2.55, -2.93, and —2.88, respectively. We caution, however, 
that the exact slope .v depends intimately upon our adopted velocity evolu- 
tion, estimated radio size, and adopted distance. The early upper Umit (red) 
is the required density jump needed to reconcile the average estimated initial 
shock velocity with the much smaller average velocities inferred from the op- 
tical spectrum, assuming the density represented by the last measured point 
is fully ionized. The dotted lines denote constant values of M/vw in units of 
Mq yr"' km"' s and bracket the typical values for red supergiant (10^^) and 
Wolf-Rayet (10^^) winds. For comparison, we also show the crude density 
profile (thin grey curves) for SN 1987A: the inner segment (left) is the density 
structure required to explain both the X-ray and radio emission (Dwarkadas 
2007c), while the outer Hnes (right) denote the upper and lower bounds on 
the density structure as deduced by light echoes (Sugerman et al. 2005). 

may indicate that a large portion of the CSM being overtaken 
is actually neutral, and thus not traced by the radio absorption. 
Since recombination is a strong function of density, this dis- 
connect may further imply that there is significant clumping 
within the CSM of SN 1996cr. 

Although we are unable to constrain the density profile en- 
countered by the blast wave at early times via Tff, the early 
implied velocity drop argues for a density increase of >120, 
while the radio luminosity jump supports a similar increase 
(even larger if oc Px^u-(^qIm holds). Thus it appears that 
we are seeing the transition between two distinctly different 
density regions separated by a dense shell, possibly associated 
with a wind-blown bubble (e.g., Dwarkadas 2005). If we as- 
sume the SN exploded into a sparse stellar wind (pcsM oc r"^) 
initially, then we might expect the shock to adhere to a self- 
similar solution prior to 1996 August (i.e., radio emission 
evolves as Ly oc f-"5-r' ^ CF94, ChevaHer 1996). Thus 
for explosion dates earlier than 1995-06-07 (i.e., day 98 for 
fo=1995-02-28), our 8.5 GHz upper limit would constrain the 
CSM density in this region to be even lower still if no fur- 
ther SSA or FF absorption is present. A lower density, how- 
ever, would imply a substantially larger density jump later on, 
which is difficult to reconcile with our current shock velocity 

unperturbed. 



evolution. This suggests that either the SN went off after this 
date, the early density did not follow a profile, or there was 
additional heavy absorption that affected even the 8.5 GHz 
data. 

At the assumed radii, the CSM density of the shells equates 
to 0.05 M© located between (1.7-2.8) x 10'^ cm. Addition- 
ally, the SN blast wave took only « 1-2 yr to reach this high- 
density region, equating to '^(l-2)vs/vw yr in the progenitor 
time frame. For reasonable wind velocities, this implies the 
region was formed a mere ^^70-14,000 yr prior to the explo- 
sion. 

For comparison, the density for a constant progenitor wind 
(i.e., constant M/vw) can be written as 

PCSM = 6.3 xl0-'2( 



M 



M. 



o yr 



km s 
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(7) 

and a few lines of constant M/v^ are plotted in Fig. 12. No- 
tably, the evolution of SN 1996cr beyond ^10'^ cm appears to 
be consistent with a constant wind of M/v„ ^ 10~^-10~'^, typ- 
ical of a Wolf-Rayet (WR) stage. However, it seems to have 
transitioned to a much sparser wind of M/v„ < lO"*^ before 
this (i.e., during its last many decades/centuries). 

There are very few scenarios that can lead to such a CSM 
density distribution, with an initial low density, followed by 
a large density jump, and then a region of decreasing density. 
Such a distribution is reminiscent, however, of a wind-blown 
bubble, wherein the interaction of a star's stellar wind with 
either a previous stage of the stellar wind or with a constant- 
density medium such as the ISM, results in material being 
swept up into a thin dense shell between the two regions. Go- 
ing outwards in radius, the bubble density distribution is de- 
scribed by (Weaver et al. 1977; Dwarkadas 2005): A region 
of freely expanding wind decreasing as r"^, the wind termi- 
nation shock, a constant-density region of shocked wind, the 
thin dense shell of swept-up ambient material and the outer 
medium, which could be another wind region with density 
decreasing as r"^. The circumstellar distribution derived from 
the radio observations seems to fit this well, with the high- 
density jump denoting the location of the dense shell, which 
we may not have seen from free-free absorption if the dense 
shell was not highly ionized. Interior to the shell is the low 
density shocked and unshocked wind, for which we only have 
upper limits. Exterior to the shell is the unshocked ambient 
medium, presumably also a wind, which we verify from our 
free-free absorption. If true, this picture again indicates that 
the stellar parameters changed significantly just a few years 
to a few thousand years before the death of the star, and that 
this region was formed by wind-wind interaction. Although 
the shell is thin, it is extremely dense, and the shock colliding 
with this region could be slowed dramatically, and lead to an 
increase in radio and X-ray emission, as is observed. 

SN 1996cr can be compared to other observed young RSNe, 
among which SN 1987A (see McCray 2005, 2007, for re- 
views) is again a good starting point because of its rising ra- 
dio light curve. SN 1987 A was detected at radio wavelengths 
when it first went off, but this initial outburst was very short- 
lived (Turtle et al. 1987) and attributed to shock acceleration 
of synchrotron-emitting electrons in the stellar wind close 
to the star at the time of the explosion (Storey & Manch- 
ester 1987; Chevalier & Fransson 1987). Radio emission 
was again detected from the supernova after ~3 yr in quies- 
cence, at first rising in a dramatic fashion for several hundred 
days and increasing more or less monotonically thereafter (see 
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Fig. 6). Unlike the early absorption in SN 1996cr, however, 
the observed spectral index of SN 1987 A has remained rela- 
tively constant between -1.0 to -0.9, indicating unabsorbed 
synchrotron emission at all times during this re-emergence 
(Manchester et al. 2002). The steady rise in SN 1987A is also 
much steeper than in SN 1996cr and > lO'* times fainter, in- 
dicating a markedly different CSM-density distribution, and 
perhaps sharper density contrasts. The constant monitoring 
at radio, optical, and X-ray wavelengths has illuminated the 
density structure surrounding SN 1987A and directly con- 
firmed multiple evolutionary stages of mass loss (Cheva- 
lier & Dwarkadas 1995; Sugerman et al. 2005; Dwarkadas 
2007a,b). Unfortunately such consistent monitoring is lack- 
ing for SN 1996cr, leaving us with only a crude understanding 
of its density structure to date. Clearly the CSM distribution 
around SN 1996cr appears to be much denser and more com- 
pact than SN 1987 A, in addition to potentially showing a pile- 
up of material around ^ 10'^ cm. We note that more modest 
structure in the CSM density, perhaps from the wind-wind in- 
teraction of multiple evolutionary stages or inhomogeneities 
within a particular stage, has been observed in many well- 
sampled RSNe light curves now, including SN 1978K (Smith 
et al. 2007a), SN 1979C (Montes et al. 2000), SN 1980K 
(Montes et al. 1998), and SN 1993J (Bartel et al. 2002; Weiler 
et al. 2007), and demonstrates that perturbations from pcsM oc 
are probably commonplace. With finer sampling in the 
future, we should be able to place stronger constraints on the 
outlying density structure in SN 1996cr. 

4.4. X-ray Light Curve 

We turn next to the X-ray light curve of SN 1996cr. The 
Circinus Galaxy was observed with a variety of X-ray instru- 
ments as shown in Fig. 6 and Table 5. From these, we find 
that SN 1996cr exhibits strong, atypical temporal evolution. 
Unlike most SNe, SN 1996cr is not detected at early times 
despite relatively frequent and sensitive observations. Since 
thermal X-ray emission effectively traces the density of the 
CSM, these early X-ray upper limits, especially when com- 
pared to the detections later on, again indicate the presence 
of either strong early absorption (A^h ^ 2 x lO'^^ cm~^) or a 
low-density cavity immediately surrounding the progenitor. 
Our radio constraints suggest in fact that both are present and 
likely to affect the X-ray emission: the mass associated with 
the FF absorption is already marginally enough to obscure 
an extrapolation of the detected X-ray slope, while the lower 
density implied by the radio luminosity jump would also re- 
sult in decreased early X-ray emission. When finally detected 
on 2000-01-16, the 0.5-2 and 2-10 keV X-ray fluxes are best- 
fitted as txt^-^'^^^ '^'' andf°^''=^'^-''', respectively. The slopes are 
consistent with each other to within the errors, suggesting that 
any apparent emerging soft component is marginal and that 
the detected portion of the rise is likely to be optically thin. 
This strong, continued rise is in stark contrast to the evolution 
of nearly all other SNe detected to date, which are theoreti- 
cally expected to decline as oc t'^-t"'^ '^ (e.g., CF94) and ob- 
servationally follow suit albeit with more scatter (e.g., see Fig. 
2 of Immler & Kuntz 2005). Sustained increases in the X-ray 
have only been observed for SN 1987 A (Park et al. 2006), 
SN 2006jc (Immler et al. 2008), and perhaps marginally for 
SN 1978K (based upon a single early upper Umit; Schlegel 
et al. 2004). 

The X-ray luminosity of the forward and reverse shocks at 
a particular radius r is given by 



Lx = 4Tr J As(Te,i)nlirfdr, (8) 

where the index ; refers to values for either the forward or 
reverse shock and Afj is the temperature-dependent cooling 
rate (Chevalier & Fransson 2003). The detection of strong 
He-like Fe and Si emission lines in the X-ray spectra (Bauer 
et al. 2001) suggests the presence of heavy elements in the 
X-ray emitting material, and therefore that the X-ray emis- 
sion stems from the shocked ejecta material found behind the 
reverse shock. 

Because the progenitor CSM distribution likely includes the 
presence of pre-existing shocks and density discontinuities as- 
sociated with the wind-blown bubble (see §4.3), the interac- 
tion of the SN blast wave with this medium is complex. Due 
to the presence of various discontinuities, the standard model 
(e.g ChevaUer 1982b, ; CF94) of the interaction of a SN blast 
wave with a power-law ambient density medium is not ap- 
plicable, and we must rely on detailed hydrodynamical mod- 
els to infer the specific mass-loss properties of the progeni- 
tor (e.g., Dwarkadas 2005, 2007a,c). Furthermore the CSM 
may be at least partially photoionized, as is thought to be the 
case for SN 1987A (e.g.. Chevalier & Dwarkadas 1995). As 
shown in the previous section for a wind-blown bubble, even 
if the progenitor star had a constant wind mass-loss rate and 
velocity (which is not likely), the structure of the CSM can 
vary substantially with radius. Thus we cannot assume that 
the CSM density directly translates into the wind parameters. 
Nonetheless, we can make some quaUtative arguments based 
on Equation 8. 

Between days 1800-4400, we find that the intrinsic X-ray 
luminosity scales almost linearly with time. The radius over 
this time period should increase by a factor of «1.6, depend- 
ing on the exact nature of the ejecta density profile and the 
CSM. The width of the emitting region can be approximated 
to be about 0.1 times the radius (e.g., CF94), resulting in 
an overall volume change by a factor of <5. From the X- 
ray spectrum, we do not see any significant changes in the 
temperature or absorption, indicating that A//(7]) remains ap- 
proximately constant. If the X-ray emission was arising in the 
shocked ambient medium, the density in Equation 8 must de- 
crease by a factor of «2.2, or equivalently be oc r"' to produce 
the observed dependence. This is significantly different from 
either a constant density or the « r~^ dependence expected 
for a stellar wind with constant parameters. Furthermore, for 
a density profile decreasing as r~' the velocity would be ex- 
pected to decrease with time, and therefore the temperature 
(which is proportional to v^) would drop accordingly. The fact 
that this is not seen makes it unlikely that the X-ray emission 
is arising from the shocked wind (see Fig. 7). 

We instead contend that the emission is arising from 
shocked ejecta behind (in a Lagrangian sense) the reverse 
shock. This scenario is strengthened by the presence of heavy 
element emission hues in the X-ray spectrum (cf. Fig. 6 of 
Bauer et al. 2001), since such elements are far more likely 
to reside in the processed ejecta than in the progenitor's stel- 
lar wind. If the SN ejecta did collide with a dense shell of 
CSM from a pre-existing wind-blown bubble, then the for- 
ward shock would advance quite slowly, and the bulk of the 
X-ray emission would arise from the reverse- shocked mate- 
rial. The post-shock temperature just behind the reverse shock 
could be quite high (>10^ K) if the density were low — 
technically the density could drop to zero for expansion in 
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a constant-density medium. The electron temperature, how- 
ever, might be an order of magnitude lower if the electrons and 
ions have not yet reached equihbrium. The temperature would 
not change appreciably over this period, consistent with our 
observational constraints. The velocity of «3400 km s~' in- 
ferred from the X-ray temperature in this case would refer to 
the electron temperature in the reverse shock, and thus the ac- 
tual shock velocity, as traced by the ion temperature, could 
be considerably higher. Note that a similar linearly increas- 
ing light curve profile was found for the hard X-ray emission 
from SN 1987A (Park et al. 2005, 2006). Dwarkadas (2007a) 
has carried out numerical simulations that suggest that this 
emission arises from the reverse shocked ejecta, as had been 
postulated by Park et al. (2006). One puzzle which remains, 
though, is how the strong the Ha and Oxygen emission-fine 
luminosities are powered. This discrepancy could be recon- 
ciled if there were holes in the CDS, such that a small fraction 
of the X-ray emission could escape while the rest went into 
exciting the strong optical line emission. 

We can again gain some further insight by comparing 
SN 1996cr to other SNe. As stated above, SN 1987A has the 
most similarities. Like SN 1996cr, SN 1987A was undetected 
at X-ray wavelengths when it first went off, thereby excluding 
an explosion directly into a slow, dense wind that is character- 
istic of a red supergiant. Notably, SN 1987A was detected 130 
days after explosion (Dotani et al. 1987; Sunyaev et al. 1987), 
although the X-ray emission here was attributed solely to ra- 
dioactivity and appeared to rise and fall over the following 
year. After this brief period, the X-ray emission went unde- 
tected for another » 3 yrs before resurfacing above the noise 
(Hasinger et al. 1996). Since this time, the X-ray emission 
has continued to rise as shown in Fig. 6, with its light curve 
well-modeled by a broken power law (Park et al. 2006). The 
reappearance of the emission at about 3 years was postulated 
as due to the interaction of the SN blast wave with a region 
of ionized wind material from the progenitor (ChevaUer & 
Dwarkadas 1995). They suggested that this region is mate- 
rial photoionized from the dense circumstellar shell that sur- 
rounds SN 1987A, formed as a result of the interaction of the 
progenitor blue supergiant wind with a pre-existing red super- 
giant wind. Although the evolution of the star may be differ- 
ent, a physical scenario broadly like that of SN 1987 A, with a 
region formed by mass-loss from a progenitor star, certainly 
seems plausible for SN 1996cr, although the contrasting light- 
curve slopes, luminosities, and radial distances clearly high- 
fight the different CSM density structure. The only other SN 
to exhibit a rise at X-ray wavelengths is SN 2006aj, a pecu- 
liar type lb. The light curve of SN 2006aj has been ascribed to 
the interaction of the blast wave with a dense shell of material 
left over from an LBV-fike outburst of the SN progenitor ~2 
years prior to the explosion (Immler et al. 2008). 

The X-ray spectrum for SN 1996cr, moreover, is modeled 
by a hot thermal plasma (kT =13.4 keV) similar to the early 
phase of SN 1993 J (Zimmermann & Aschenbach 2003) but 
substantiaUy hotter than SN 1987A (kT « 3 keV with an 
emerging 0.3 keV component), and implies that the condi- 
tions within the X-ray-emitting region are somewhat different. 
The hotter temperature could indicate a higher shock veloc- 
ity, or a larger level of equilibration between the electrons and 
ions. The X-ray spectrum of SN 1996cr is additionally ab- 
sorbed below ~2 keV. While some fraction of the substantial 
non-Galactic column density could be due to self-absorption, 
we consider this to be minimal given that other X-ray point 
sources in the Circinus Galaxy have similar X-ray derived col- 



umn densities (Bauer et al. 2001). 

The rapid evolution we see here for SN 1996cr has been 
noted in a few other objects, such as SN 2002kg (Van Dyk 
et al. 2006; Maund et al. 2006) and SN2006jc (Foley et al. 
2007; Pastorello et al. 2007). Such an evolution has also 
been invoked to explain the enormous energy output from 
SN2006gy (Ofek et al. 2007; Smith et al. 2007b), although 
the lack of obvious CSM-interaction measures at X-ray, Ha, 
and radio wavelengths here implies that perhaps something 
more extreme is going on (e.g., Smith & McCray 2007c). 
Unfortunately, for many of the older, classic X-ray emit- 
ting type lln SNe such as SN 1978K (Schlegel et al. 2004; 
Smith et al. 2007a), SN 1979C (Immler & Kuntz 2005), and 
SN 1986J (Temple et al. 2005), early-time behavior of the X- 
ray emission remains completely unknown and thus the fre- 
quency and duration of sudden evolutionary changes in SNe 
progenitors is still an open question. The X-ray evolution of 
SN 1996cr may allow us to explain the relatively flat late-time 
X-ray ligfit curves for objects fike SN 1978K and SN 1979C, 
for instance. 

5. CONCLUSIONS 

We have confirmed SN 1996cr as one of the nearest, and 
subsequently X-ray and radio-brightest, SNe to date. Archival 
optical imaging data have constrained the explosion date to 
within Rsl yr prior to 1996-03-16, while a new VLT spec- 
trum has identified it as a type lln SN. The multi-wavelength 
constraints detailed above point to a rich progenitor history. 
As is the case for most type lln SNe, SN 1996cr appears to 
have originated from a massive star that shed its outer layers 
at some late evolutionary stage, which now act as the target 
for the outgoing SN shock. 

The optical light curve hints at a mild re-brightening at late 
times, possibly driven by line emission, or by X-ray photons 
being downscattered and being emitted as optical photons. 
From the recent optical spectrum, the strong narrow H emis- 
sion points to a significant accumulation of mass in the CSM 
prior to the SN. We are additionally able to resolve the com- 
plex O emission into a number of distinct high-velocity com- 
ponents wifli FWHM~2000-3000 km s"^ The origin of such 
distinct O complexes remains somewhat unclear, perhaps aris- 
ing from one or more of the following: structural asymme- 
tries, differential absorption, or multiple concentric shocks. 

The X-ray and radio emission allow us to quantify the CSM 
density, tracing out two apparently distinct regions, which 
cannot be reconciled with a single constant-wind model. 
From the radio and X-ray upper limits, we infer that the in- 
ner region (r < 1 x 10'^ cm) is likely quite sparse. The outer 
region [r « (1.5-2.8) x lO'^ cm], as traced by the early radio 
absorption, radio luminosity jump, and X-ray emission, is not 
well constrained, and is Ukely associated with the interaction 
region between two distinct progenitor winds, which culmi- 
nates in a thin, dense wind-swept shell. This second region 
appears relatively dense (pcsM ~ 10^ cm"^), potentially con- 
sistent with a slow- velocity WR-fike wind, although the CSM 
density profile traced by tiie radio absorption (pcsM oc r~^^) 
naively seems at odds with the implied CSM density pro- 
files from the radio and X-ray luminosities (pcsM occonstant 
or r"', respectively). A more detailed understanding of the 
CSM awaits sophisticated hydrodynamical models. 

The most plausible explanation is one where the progeni- 
tor changed evolutionary states just prior to explosion, likely 
transitioning from a WR-like wind to an even faster, less 
dense wind, resulting in portions of the previous CSM be- 
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ing swept-up to form a wind-blown cavity surrounded by a 
dense shell. Such a scenario has been recently proposed for 
the peculiar type lb SN2006aj (Foley et al. 2007; Pastorello 
at al. 2007; Immler et al. 2008). The radio and X-ray emis- 
sion are constant or still rising, demonstrating that the blast 
wave is still making its way through the dense, swept-up shell 
at the edge of the cavity. The collision of the blast wave with 
the shell would result in a transmitted shock expanding into 
the shell, and a reflected shock expanding back into the SN 
ejecta. This reflected shock should be distinguished initially 
from the reverse shock of the SN expanding into the ejecta, 
although it will eventually overtake the reverse shock. The 
overall effect could be one of several shocks and rarefaction 
waves travelling back and forth through the ejecta. A spec- 
trum would reveal several different velocity structures, per- 
haps as we observe in the O lines. The reason why it is seen 
mainly in O is unclear, but it could indicate that this is an O- 
rich remnant, similar to Cas A for example (see Fesen et al. 
2001, and references therein). Or perhaps it is just that the 
reverse shock has not penetrated very far into the ejecta, so 
that only the O-rich layers have been shocked. We note that if 
the shock is indeed progressing into shell of constant density, 
then the resulting density gradients would not lead to the for- 
mation of a CDS at all (e.g., see Fig. 1-4- of Chevalier 1982a). 
If so, then the observed O velocity components may instead 
arise from a disk or ring prehaps similar to SN 1987A. 

Once the blast wave makes its way through the dense shell, 
it will emerge from the other side, presumably as a weak 
shock. This should encounter a relatively pristine wind and 
start to decUne in both the radio and X-ray bands. Further 
regular monitoring of SN 1996cr should reveal this. Once this 
happens, we will be able to strengthen many of our current 
constraints and assumptions. 

We have compared SN 1996cr to the late-time behavior of 
SN 1987A, which is the only other SN with a documented 
dramatic increase at X-ray and radio wavelengths. SN 1996cr 
is likely propagating through its CSM in a manner similar 
to SN 1987 A expanding into its H II region and circumstel- 
lar ring (e.g.. Park et al. 2006). The X-ray and radio lumi- 
nosities of SN 1996cr, however, are several orders of magni- 
tude brighter and comparable to typical type Iln SNe such 
as SN 1978K and SN 1979C. This is presumably because 
the swept-up CSM is more compact, and perhaps, unlike 
SN 1987 A, there is no intervening high-density material en- 
countered which slows down the ejecta until they reach a 
dense shell, so the density jump is higher and subsequently 
more dramatic. As we follow SN 1996cr over the next few 
decades, we anticipate that it will transition into an evolution 
similar to other SNe. It is intriguing to think that SN 1996cr 
might represent a sort of "bridge" object between the extremes 
of SN 1987A and these more luminous type II SNe. Although 
plagued by small number statistics, the fact that two (al- 
though perhaps more if we count SN 1978K and SN 1979C) 
of the ^5-10 closest SNe to have exploded in the past four 
decades show evidence for wind-blown bubbles impUes that 
the phenomenon is conmion, and may indicate that current 
SNe searches are somehow biased against them at larger dis- 
tances. Systematic constraints exist only for type Ibc SNe 
(Soderberg et al. 2006a,b), where late-time radio observations 
indicate that only ~2% of initially undetected targets show 
late-time detections, while ~40% of detections exhibit abrupt 
light-curve variations. 

The large mass associated with the pre-SN CSM sur- 
rounding SN 1996cr could have been the result of explosive 



mass loss thousands of years before the SN explosion (e.g., 
Woosley et al. 2002; Heger & Woosley 2002), qualitatively 
similar to SN 2006aj (Foley et al. 2007; Pastorello et al. 2007; 
Immler et al. 2008), although the limited X-ray and radio data 
point to a relatively smooth density profile, which implies a 
more gradual process was at work. Numerous other SNe in- 
cluding SN 1994W (Chugai et al. 2004), SN2001em (Chugai 
& Chevalier 2006; Bietenholz & Bartel 2007), SN2006gy 
(Ofek et al. 2007; Smith et al. 2007b), and SN2006jc (Fo- 
ley et al. 2007; Pastorello et al. 2007; Inomler et al. 2008) all 
show signs of dense, enriched CSM likely produced by gi- 
ant mass-loss events just prior to their SNe, implying that this 
form of mass-loss is relatively common. LBVs have often 
been invoked to account for such episodes of extreme mass- 
loss (e.g., Gal- Yam et al. 2007; Smith et al. 2007b). Indeed a 
giant outburst was seen in the case of SN 2006jc just 2 years 
prior to core-collapse (Pastorello et al. 2007). Again it is un- 
clear whether the explosive nature of LBVs can be recon- 
ciled with the relatively smooth mass-loss/density gradients 
observed for SN 1996cr. 

Our analyses clearly demonstrate that SN1996cr is a 
unique and compelling target. Its proximity and flux enable 
a wealth of observations that would be otherwise impracti- 
cal for typical SNe. We have reported on initial findings of 
a VLBI campaign underway to resolve spatially SN 1996cr 
at several radio wavelengths, eventually enabling studies of 
its morphology and temporal expansion. This will provide a 
more coherent picture, allowing us to set strong limits on the 
expansion velocity from the resolved size and potentially con- 
firm the impUed asymmetries we infer from the optical spec- 
trum. Follow-up HST imaging will provide constraints on any 
potential light echoes in the vicinity of SN 1996cr, possibly 
allowing refined estimates of its explosion date, distance, and 
outlying CSM structure. Likewise, follow-up X-ray observa- 
tions will hopefully elucidate the nature of the strong contin- 
ued rise in the X-ray band and better constrain properties of 
the reverse shock region. Notably, SN 1996cr is bright enough 
to be efficiently monitored using the Chandra HETGS, en- 
abling high-resolution studies of its numerous X-ray emission 
lines and their potential evolution; an in-depth discussion of 
the current X-ray spectral constraints, including a compos- 
ite Chandra HETGS spectrum from the data presented here, 
are forthcoming (F. Bauer et al., in preparation). Addition- 
ally, future optical spectroscopic studies should allow investi- 
gations into the temporal behavior of complex emission lines 
found in our discovery spectrum. Infrared imaging and spec- 
troscopy will help address the potential for dust formation; 
notably type 11 SNe are thought to contribute strongly to the 
overall dust content of galaxies. Looking ahead to the near 
future, if SN 1996cr remains relatively bright at milhmeter 
wavelengths, we anticipate that observations with ALMA will 
provide a substantial leap in terms of understanding the struc- 
ture and composition of SN 1996cr and its relation to both 
SN 1987A and more typical type 11 SNe. 

If any researchers have further serendipitous multi- 
wavelength observations of SN1996cr during the critical 
1995-1996 period, we would be grateful to learn of these. 
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